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ABSTRACT 


A  technique  for  the  use  of  a  digital  computer  for  simul- 
taneous generation  of  feedback  control  system  transient 
response  and  sensitivity  to  parameter  variations  is  developed. 
The  sensitivity  model  for  a  basic  linear  control  system  is  a 
cascaded  replica  of  the  original  system  with  the  sensitivity 
pickoff  points  corresponding  to  the  error  and  feedback  signals 
of  the  original  system.   This  is  convenient  for  computer  anal- 
ysis with  DSL/360.   For  systems  with  nonlinear  elements,  the 
sensitivity  model  is  implemented  by  making  the  instantaneous 
gain  equal  to  the  slope  of  the  nonlinear  function.   Sensi- 
tivity analysis  is  utilized  for  improvement  of  system 
performance,  determining  parameter  tolerances,  and  predicting 
changes  in  the  response  due  to  parameter  variations. 


Library 

U.S.    Naval  Postgraduate  School 

Monterey    California       93940 


TABLE    OF    CONTENTS 

Page 

I.   INTRODUCTION  9 

II.   BASIC  CONTROL  SYSTEM  SENSITIVITY  12 

A.  INTRODUCTION 12 

B.  GAIN-CONSTANT  SENSITIVITY  COEFFICIENTS  13 

C.  APPLICATION i6 

III.   MULTILOOP  CONTROL  SYSTEM  SENSITIVITY 25 

A.  INTRODUCTION 25 

B.  GAIN-CONSTANT  SENSITIVITY  COEFFICIENTS  25 

C.  TIME-CONSTANT  SENSITIVITY  COEFFICIENTS 30 

D.  APPLICATION 31 

IV.   EXAMPLES  OF  DIFFERENTIAL  EQUATION  SENSITIVITY 

ANALYSIS 41 

A.  INTRODUCTION 41 

B.  SECOND-ORDER  DIFFERENTIAL  EQUATION 
SENSITIVITY  TO  COEFFICIENT  VARIATIONS  41 

C.  SIXTH-ORDER  DIFFERENTIAL  EQUATION  SENSITIVITY 
TO  COEFFICIENT  VARIATIONS  43 

V.   NONLINEAR  SYSTEM  SENSITIVITY  53 

A.  INTRODUCTION 53 

B.  SENSITIVITY  COEFFICIENT  FOR  A  GAIN  FACTOR  5  3 

C.  APPLICATION 56 

VI.   CONCLUSIONS 65 

APPENDIX  A  -  DSL/360  DIGITAL  SIMULATION  PROGRAM  FOR 

BASIC  SYSTEM 66 

APPENDIX  B  -  DSL/360  DIGITAL  SIMULATION  PROGRAM  FOR 

MULTILOOP  SYSTEM  67 


TABLE  OF  CONTENTS  (Continued) 

Page 

APPENDIX  C  -  DSL/360  DIGITAL  SIMULATION  PROGRAM  FOR 

SIXTH-ORDER  DIFFERENTIAL  EQUATION  68 

APPENDIX  D  -  DSL/360  DIGITAL  SIMULATION  PROGRAM  FOR 

NONLINEAR  SYSTEM  69 

LIST  OF  REFERENCES 70 

INITIAL  DISTRIBUTION  LIST 71 

FORM  DD  1473  73 


LIST  OF  FIGURES 

FIGURE  Page 

2.1  Canonical  form  of  a  feedback  control  system 14 

2.2  Feedback  control  system  and  sensitivity  model 14 

2.3  Basic  control  system  response  and  sensitivity 

for  K  =  50  15 

2.4  Overshoot  reduction  of  0.1  by  utilization  of 
sensitivity  coefficients 19 

2.5  Gain  tolerance  of  14%  for  overshoot  of  +  0.05  ----20 

3.1  Multiloop  system 27 

3.2  Multiloop  sensitivity  model 29 

3.3  Multiloop  system  response  and  sensitivity 
functions 32 

3.4  Amplifier-exciter  sensitivity  functions 33 

3.5  Multiloop  response  for  K  =  40  and  K  =  60 33 

3.6  Load-  and  feedback-gain  sensitivity  functions 34 

3.7  Multiloop  response  for  K  =  1.1  and  F  =  2.2 34 

4.1  Block  diagram  of  mass-spring-damper  system 42 

4.2  Second-order  differential  equation  sensitivity 
model 42 

4.3  Sixth-order  differential  equation  sensitivity 

model 44 

4.4  Second-order  system  response  and  sensitivity 
functions 45 

4.5  Sixth-order  unstable  system  response  and 
sensitivity  functions  46 

4.6  Sixth-order  stable  response  and  sensitivity 
functions 48 

5.1    Block  diagram  of  nonlinear  system  54 


LIST  OF  FIGURES  (Continued) 

FIGURE  Page 

5.2  Sensitivity  model  of  nonlinear  system 54 

5.3  Characteristics  of  nonlinear  element  N 58 

5.4  Characteristics  of  sensitivity  model 
nonlinearity 58 

5.5  Response  and  sensitivity  for  nonlinear  system 
with  pi  =  -0.1  and  p2  =  0.1 59 

5.6  Response  and  sensitivity  for  unlimited  linear 
system 59 

5.7  Response  and  sensitivity  for  nonlinear  system 
with  pi  =  -0.2  and  p2  =  0.2 61 

5.8  Gain  tolerance  of  14%  for  overshoot  of  +  0.014-  62 


LIST  OF  TABLES 

TABLE  Page 

2.1  Computer  output  data  for  unity  feedback  system  21 

2.2  Computer  output  data  for  K  =  55,  50,  45,  and 

36 23 

3.1a   Computer  output  data  for  C,  S„  ,  and  S_   37 

KL  1  i 

3.1b    Computer  output  data  for  C,  S„  ,  S_ ,  and  S„  39 

K3    F        TH 

4.1    Computer  output  data  for  X,  S_ ,  and  S,, 51 

5.1  Comparison  of  actual  outputs  and  outputs 

predicted  from  sensitivity  coefficients  61 

5.2  Computer  output  data  for  nonlinear  system 

outputs  and  sensitivity  63 


I.   INTRODUCTION 

The  early  development  of  sensitivity  analysis  was 
applied  to  circuit  theory  where  it  is  useful  to  know  how 
sensitive  the  response  of  the  circuit  is  to  variations  of 
circuit  parameters  such  as  resistance,  capacitance,  and 
inductance.   There  are  many  different  definitions  of  sensi- 
tivity available.   H.  Bode  [1]  defines  a  logarithmic  or 

normalized  sensitivity  function  of  a  variable,  x. (t,a.) , 

x.    Hn  a. 
with  respect  to  a  parameter,  a .  ,  as  S    =  Sj,  '  •   ^    and  I.  M. 

3        j      n   i 

Hortowitz  [2]  adapted  the  inverse  of  this  definition  which 
is  widely  used.   J.  Leeds  [3]  presented  sensitivity  function 
for  incremental  changes  in  parameters  and  S.  R.  Parker  [4] 
recently  generalized  the  technique  of  sensitivity  analysis  for 
nonlinear  circuits.   Much  of  this  work  entailed  the  definition 

of  R.  Tomovic  [5]  who  defined  the  sensitivity  coefficients  as 

x.    3x.  x.    8x. 

S„   =  * .  Kokotovic  [6]  adapted  the  definition  S    =  yss 

a. 3a.  r  a.d£na. 

1  3  3         3 

for  use  in  control  systems  and  developed  a  method  of  generat- 
ing sensitivity  coefficients  using  analog  computers.   In  this 
thesis  the  method  of  Kokotovic  [61  is  expanded  for  use  with 
digital  computer  techniques  utilizing  block  diagrams  and 
transfer  functions  of  forward  and  feedback  loops  in  control 
systems.   The  IBM/360  Digital  Simulation  Language  (DSL) 
program  [7]  is  very  well  suited  for  modeling  a  control  system 
and  obtaining  the  time  response  and  sensitivity  coefficients 
simultaneously . 


DSL  is  a  continuous  system  simulator  whose  non-procedural, 
problem-oriented  input  language  has  the  attractive  building- 
block  approach  of  digital-analog  simulators  with  the  added 
power  of  logical  and  algebraic  equation  notation.   The  DSL 
system  provides  a  basic  set  of  35  blocks  from  which  a  phys- 
ical system  may  be  modeled.   Included  in  this  set  are  such 
conventional  analog  components  as  integrators,  comparators, 
and  transfer  functions.   The  Kokotovic  [6]  method  of  sensi- 
tivity analysis  is  readily  adaptable  to  the  DSL  program  as 
control  systems  in  block-diagram  form  can  be  programed  in 
the  s-domain. 

Section  II  of  this  thesis  applies  the  Kokotovic  [6] 
method  of  generating  sensitivity  coefficients  of  the  basic 
feedback  control  system  for  the  forward  and  feedback  gain 
parameters.   Sensitivity  functions  are  derived  and  applied 
to  a  third-order  system.   Examples  are  given  for  the  use 
of  sensitivity  coefficients  in  analysis,  parameter  adjust- 
ment for  improvement  of  system  performance,  and  parameter 
tolerance  considerations.   Section  III  applies  the 
Kokotovic  [6]  method  of  sensitivity  analysis  to  multiloop 
systems.   Gain  and  time-constant  sensitivity  coefficients 
are  derived  and  applied  to  a  regulated  power  generator  as  an 
example.   Section  IV  demonstrates  the  Kokotovic  [6]  method 
for  the  general  case  of  differential  equations,  and  Section 
V  is  devoted  to  the  nonlinear  system. 

Sensitivity  analysis  of  a  control  system  offers  an 
analytical  method  for  determining  those  parameters  of  a 
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dynamic  system  which  have  the  greatest  influence  on  the 
system  performance.   High-speed  digital  computers  give 
very  precise  solutions  for  the  system  equations,  but  it  is 
impossible  to  realize  exact  mathematical  models  and  para- 
meter values.   Additional  variations  in  parameters  may  also 
be  caused  by  wear,  aging,  shifting  of  operating  point, 
temperature,  and  etc.   It  is  beneficial  to  know  the  effect 
of  these  small  parameter  variations  on  system  performance. 
Sensitivity  analysis  can  give  this  necessary  insight  so 
that  all  parameters  of  a  dynamic  system  are  not  given  equal 
importance.   This  is  more  efficient  than  repetitive  computer 
solutions  with  parameters  changed  by  increments.   Analysis 
utilizing  sensitivity  coefficients  can  put  emphasis  on 
specific  component  tolerances  and  points  that  are  critical 
to  the  performance  specifications. 
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II.   BASIC  CONTROL  SYSTEM  SENSITIVITY 

A.   INTRODUCTION 

The  Kokotovic  [6]  method  for  computing  sensitivity 
coefficients  for  various  parameters  is  restricted  to  linear, 
time-invariant,  deterministic  systems.   Although  it  is 
true  that  physical  systems  do  not  have  these  ideal  properties, 
linear  approximations  for  these  systems  have  proven  adequate 
in  many  useful  cases.   When  parameter  values  are  to  be 
approximated,  it  is  important  to  estimate  the  influence  of 
each  parameter  on  the  system  response.   This  can  be  done 
accurately  by  sensitivity  analysis.   Also,  limits  of  indi- 
vidual parameter  variations  may  be  deduced  and  tolerance 
levels  formulated. 

The  Kokotovic  [6]  method  of  solving  for  sensitivity 
coefficients  has  the  advantage  that  all  of  the  functions 
are  generated  simultaneously  with  the  time-domain  solution. 
Kokotovic  [6]  gives  analytic  proof  that  this  is  possible. 
As  an  engineering  approach  to  analysis,  it  is  common  to 
approximate  a  high-order  system  by  an  equivalent  second- 
order  system  to  enable  using  graphs  and  tables  for  second- 
order  response  characteristics.   Qualitative  observations 
of  graphs  of  sensitivity  coefficients  can  give  insight  into 
the  transient  performance  in  terms  of  overshoot,  settling 
time,  rise  time,  and  etc.,  for  higher  degree  characteristic 
equations  as  will  be  seen  in  examples, 
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B.   GAIN  CONSTANT  SENSITIVITY  COEFFICIENTS 

The  sensitivity  coefficient  S„  is  defined  as  the  ratio 
of  the  change  in  output  to  the  fractional  change  in  the 
parameter  K. 

S   =  ?C  '   =   9C  (7    1) 

bK  ~  W*      9£nK  ^'X) 

where  C  is  the  controlled  output  and  K  is  the  parameter  being 

investigated. 

The  canonical  form  of  a  feedback  control  system  after 

Laplace  transformation  is  shown  in  Fig.  2.1.   The  definitions 

of  terms  are  as  follows: 

K  =  direct  gain  constant 

KG(s)  =  direct  transfer  function 

F  =  feedback  gain  constant 

FH(s)  =  feedback  transfer  function 

KG(s)FH(s)  =  loop  transfer  function 

R(s)  =  reference  input 

E  =  error  signal 

B  =  feedback  signal 

The  basic  equations  are  as  follows: 

R{    '  l+KG(s)FH(s)  U.^aj 

R(S)   =   l+KG(s)FH(s)  (2.2b) 

B         KG(s)FH(s) 

RKS)  l+KG(s)FH(s)  U.ZC) 

In  the  work  that  follows,  capital  letters  denote  the 
Laplace  transformation  of  the  functions  described  except  for 
the  letters  K  and  F  which  are  constants. 
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Fig. 2.1.  Canonical  form  of  a  feedback  control  system 
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Fig. 2. 2.  Feedback  control  system  and  sensitivity  model. 
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Fig. 2. 3.  Basic  control  system  response  and  sensitivity 
for  K  =  50. 
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Applying  (2.1)  to  (2.2a)  yields  the  gain  sensitivity 

coefficients 

RKG  n 


SK   " 

(1+KGFH) 2 

1+KGFH 

SF   = 

-RK2G?FH 
(1+KGFH) 2 

C(-KGFH) 
1+KGFH 

(2.3) 


(2.4) 


These  two  gain  sensitivity  functions  are  shown  in  block 
diagram  in  Fig,  2.2.   This  shows  that  the  sensitivity  model 
is  an  exact  replica  of  the  original  system  with  the  controlled 
output  as  the  input  to  the  sensitivity  model.   The  error  and 
feedback  signals  of  the  sensitivity  model  are  the  desired 
sensitivity  functions.   It  should  be  noted  that  the  sensi- 
tivity model  requires  no  new  components  or  transfer  functions 
and  all  quantities  axe  generated  simultaneously. 

C.   APPLICATION 

In  this  thesis  the  equations  for  the  feedback  systems 
are  programed  on  the  IBM  360  computer  using  the  Digital 
Simulation  Language  (DSL)  program  for  both  the  transient 
solution  and  the  sensitivity  coefficients.   This  program 
has  the  advantages  of  programing  in  the  s -domain  in  trans- 
fer function  form  and  the  accuracy  of  the  digital  computer 
with  graphical  and  numerical  results.   A  sample  program  is 
included  in  Appendix  A  and  the  numerical  solution  and 
sensitivity  coefficients  are  listed  in  Table  2.1c   Fig.  2.3 
is  the  graph  of  the  output  and  sensitivity  coefficients. 
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The  direct  transfer  function  is 

s (s+4) (s+6) 

For  example,  data  from  Table  2.1  when  S„  is  a  maximum 

Time         T   =   1.0 

Gain         K   =   50 

Sensitivity  S__  =   0.70863 
K 

Input        R   =   unit-step 

Output       C   =   0.96  45  6 

For  large  variations  (2.1)  becomes 

SK   .   «^  (2.6) 

Let  the  gain,  K,  change  10%,  then 

AC  -  S„  ~     =   0.70863(0.1)   =   0.070863 

i\   is. 

This  shows  that  the  output  will  change  approximately  +  0.071 

if  the  gain  varies  from  50  to  55  or  45.   This  can  be 

verified  from  Table  2.2  for  the  following  data: 

K  =  55 

T  =  1.0 

C  =  1.0331 

C  =  C  +  AC  =  0.96456  +  0.070863  =  1.03541 

and  for 

K  =  45 

T  =  1.0 

C  =  0.89133 

C  =  C-AC  =  0.96456  -  0.07086  =  0.89370 

This  shows  that  sensitivity  coefficients  for  changes  of  gain 

as  large  as  10%  produce  accuracy  within  0.3%  when  compared  to 

incremental  intervals  as  denoted  in  partial  derivatives. 
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Another  example  for  the  use  of  sensitivity  coefficients 
is  to  compute  the  gain  required  to  reduce  the  peak  overshoot 
of  the  system  transient  response  to  a  more  desirable  specifi- 
cation of  one-half  its  initial  value.   Table  2.1  gives  the 
following  data: 

K   =   50 

Peak  overshoot   =   0.2 

T   =   1.5 

SR  =  0.36568 

Rearranging  (2.6) 

ACK  ,_  _, 

AK   a  q (2.7) 

K'   =  K  -  AK  =  K(l-|£)  (2.8) 

bK 

=   50  (1  "  OM»  *  36 

Table  2.2  for  K  =  36  shows  that  the  peak  overshoot  is  reduced 
to  0 . 1  and  Fig.  2.4  shows  the  reduction  of  peak  overshoot  by 
a  factor  of  one-half. 

Any  other  design  specification,  such  as  settling  time  or 
rise  time,  can  be  adjusted  by  utilization  of  the  sensitivity 
coefficient.   In  the  example  it  should  be  noted  that  reducing 
the  overshoot  by  reduction  in  gain,  increases  the  rise  time 
by  0.4  second. 

Tolerance  considerations  for  parameters  are  encompassed 
within  the  definition  of  sensitivity  coefficient. 

Rearranging  (2.6) 

Tolerance   =   AK/K   -   4^  (2.9) 

bK 
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Fig. 2. 4.   Overshoot   reduction  of  0.1  by  utilization  of 
sensitivity  coefficients. 
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Pig. 2. 5.  Gain  tolerance  of  V\$   for  overshoot  of  10.05. 
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BASIC  SERVO  SENSITIVITY  COFFFICIENTS 


TIME 

0.0 

l.OOOOE-01 
2.00C0E-01 
3.0000E-01 
4.0300E-C1 
5.0000E-01 
6.0000E-01 
7.0000E-01 
8.00C0E-01 
9.0C0OE-O1 
l.OOOOE  00 
1.1000E  00 
1.2000E  00 
1.3000E  00 
1.4000E  00 
1.5000E  00 
1.600GE  00 
1.7000E  00 
1.8000E  00 
1.9000E  00 
2.0000E  00 
2.1000E  00 
2.2000E  00 
2.3OC0E  00 
2.4000E  00 
2.5000E  00 
2.6000E  00 
2.7000E  00 
2.8O00E  00 
2.900CE  00 
3.0000E  00 
3.1CC0E  00 
3.2000E  00 
3.300CE  00 
3.4000E  00 
3.5000E  00 
3.600CE  00 
3.7000E  00 
3.8000E  00 
3.9C00E  00 
4.0000E  00 
4.1000E  00 
4.2000E  00 
4.3000E  00 
4.4000E  00 
4.5000E  00 
4.6000E  00 
4.7000E  00 
4.8000E  00 
4.9000E  00 
5.00C0E  00 
5.100CE  00 
5.2000E  00 
5.3000E  00 
5.4000E  00 


OUT 
0.0 

6.5182E-03 

4.1385E-02 

1.1165E-01 

2.1292E-01 

3.3623E-31 

4.7147E-01 

6.0911E-01 

7.41C5E-C1 

8.6100E-01 

9.6456E-01 

.0491E    00 

.U37E    00 

.1586E    00 

.1853E    00 

.1959E    00 

.1931E    00 

.1797E    00 

.1587E    00 

.1328E    00 

.1046E    00 

.0760E    00 

.0489E    00 

.0245E    00 

.0036E    00 

9.8672E-01 

9.7404E-01 

9.6544E-01 

9.6056E-01 

9.5895E-01 

9.6002E-01 

9.6318E-01 

9.6782E-01 

9.7336E-01 

9.7932E-01 

9.8525E-01 

9.9083E-01 

9.9581E-01 

l.OOOOE    00 

1.0034E    00 

1.0058E    00 

1.0C75E    OC 

1.0084E    00 

1.0086E    00 

1.00  83E    00 

1.0075E    00 

1.0065E    00 

1.0053E    00 

1.0041E    00 

1.0028E    00 

1.0017E    00 

1.0007E    00 

9.9984E-01 

9.9918E-01 

9.9869E-01 


SENSK 
CO 

6.516CE-C3 
4.  1258E-02 
1.1054E-C1 
2.0812E-C1 
3.2212E-C1 
4.3885E-C1 
5.4526F-C1 
6.3042E-01 
6.864CE-C1 
7.0863E-C1 
6.959CE-01 
6.5CC5E-C1 
5.7542E-01 
4.7816E-01 
3.6558E-01 
2.4539E-01 
1.251C E-01 
1.1496E-C2 
-8.9767E-02 
-1.7441E-01 

-2.3964E-01 

-2.8418F-01 

-3.0811E-01 

-3.1273E-C1 

-3.0C32E-01 

-2.7383E-01 

-2.3663E-C1 

-1.9225E-C1 

-1.4414E-C1 

-9.5469E-C2 

-4.8996F-C2 

-6.9429E-03 

2.9C47E-02 

5.7914E-C2 

7.9165E-C2 

9.2808E-C2 

9.9281E-02 

9.9353F-C2 

9.40  37E-C2 

8.4486F-C2 

7.1906F-C2 

5.7483E-C2 

4.2310E-C2 

2.735CE-C2 

1.3390E-C2 

1.0386E-03 

-9.2916E-0  3 

-1.7368E-02 

-2.3124E-^2 

-2.6629E-02 

-2.8C75E-02 

-2.7730E-C2 

-2.5919F-C2 

-2.2996E-02 


SENSE 
-0.0 

-2.2129E-06 
-1.2685E-C4 
-1.1099E-03 
-4.7931E-03 
-1.4109E-02 
-3.2619E-02 
-6.3846E-02 
-1.1063E-P1 
-1.7460E-01 
-2.5593E-01 
-3. 5322E-C1 
-4.6363E-01 
-5.8319E-01 
-7.0712E-C1 
-8.3030E-C1 
-9.4767E-01 
-1.0546E  CO 
-1.1472E  00 
-1.2226E  00 
-1.2790E  OC 
-1.3157E  00 
-1.3331E  OC 
-1.3326E  OC 
-1.3163E  00 
-1.2870E  OC 
-1.2479E  00 
-1.2021E  00 
-1.1528E  OC 
-1.1031E  OC 
-1.0555E  OC 
-1.0122E  00 
-9.7476E-01 
-9.4432E-01 
-9.2140E-C1 
-9.0609E-01 
-8.9802E-01 
-8.9653E-01 
-9.C>66E-01 
-9.0933E-01 
-9.2136E-01 
-9.3558E-01 
-9.5088E-01 
-9.6627E-C1 
-9.8091E-CI 
-9.9414E-C1 
-1.0055E  00 
-1.0146E  OC 
-1.0214E  OC 
-1.0260E  OC 
-1.0283E  00 
-1.0288E  ^0 
-1.0276E  OC 
-1.0251E  r»C 
-1.0217E  OC 


Table   2.1.    Computer  output  data  for  unity  feedback   system 
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BASIC  SERVO  SENSITIVITY  COEFFICIENTS 


TIME 
5.5^>00E  00 
5.6000E  00 
5.7C00E  00 
5.8000E  00 
5.9000E  00 
6.0000E  00 
6.1000E  00 
6.200CE  00 
6.3000E  00 
6.4000E  00 
6.5C00E  00 
6.6000E  00 
6.7000E  CO 
6.8000E  00 
6.9C00E  00 
7.0000E  00 
7.1000E  00 
7.2000E  00 
7.3000E  00 
7.4000E  00 
7.500CE  00 
7.6000E  00 
7.7C0OE  00 
7.8O00E  00 
7.900CE  00 
8.0000E  00 
8.100CE  00 
8.2000E  00 
8.3C00E  00 
8.4OC0E  00 
8.5000E  00 
8.6O0CE  00 
8.7000E  00 
8.8000E  00 
8.9000E  00 
9.0000E  00 
9.1000E  00 
9.2000E  00 
9.3000E  00 
9.4000E  00 
9.5000E  00 
9.6000E  00 
9.7C00E  00 
9.8000E  00 
9.900CE  00 
l.OOOOE  01 


OUT 
9.9838E-01 
9.9823E-01 
9.9821E-01 
9.9829E-01 
9.9846E-01 
9.9868E-01 
9.9893E-01 
9.9920E-01 
9.9945E-01 
9.9969E-01 
9.9989E-01 
1.0001E  00 
1.0002E  00 
1.0003E  00 
1.0003E  00 
1.0C04E  00 
1.0004E  00 
1.0003E  00 
1.0003E  00 
1.0003E  00 
1.0002E  00 
L.0002E  00 
1.0001E  00 
1.0001E  00 
l.OOOOE  00 
9.9998E-01 
9.9995E-01 
9.9994E-01 
9.9992E-01 
9.9992E-01 
9.9992E-01 
9.9993E-01 
9.9993E-01 
9.9994E-01 
9.9995E-01 
9.9997E-01 
9.9998E-01 
9.9999E-01 
9.9999E-01 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 


SENSK 

1.9314E-02 

-1.52C8E-02 

■1.098CE-02 

•6.8653E-03 

-3.1281E-03 

1.4061E-04 

2.8257E-03 

4.8800E-03 

6.3002E-C3 

7.1194E-C3 

7.3986E-03 

7.2171E-03 

6.6686E-03 

5.8481E-03 

4.8504E-03 

3.7631F-03 

2.6627E-C3 

1.6125E-C3 

6.6310E-04 

-1.5C68E-C4 

-8.0872E-04 

-1.3018E-C3 

-1.6327E-C3 

-1.8120E-C3 

-1.8587E-03 

-1.7906E-C3 

•1.636CE-03 

-1.4186E-03 

-1.1616E-C3 

-8.8680E-C4 

■6.1244E-C4 

-3.5483E-C4 

1.2463E-C4 

6.9797E-C5 

2.2489E-C4 

3.3873E-C4 

4.1258E-C4 

4.4990E-C4 

4.5526E-04 

4.3404E-04 

3.9238E-04 

3.3724E-04 

2.7323E-04 

2.0635E-04 

1.4C67E-04 

8.C466F-05 


SENSF 
-1.0177E  00 
-1.0134E  00 
-1.0092E  00 
-1.0052E  or 
-1.0016E  or 
-9.9854E-C1 
-9.9611E-01 
-9.9432E-01 
-9.9315E-01 
-9.9257E-01 
-9.9249E-01 
-9.9284E-01 
-9.9352E-C1 
-9.9444E-01 
-9.9549E-01 
-9.9661E-01 
-9.9771E-01 
-9.9873E-C1 
-9.9965E-01 
-1.0004E  00 
-1.0010E  00 
-1.0015E  00 
-1.0017E  00 
-1.0019E  CC 
-1.0019E  OC 
-1.0C18E  00 
-1.0016E  Of 
-1.0014E  CO 
-1.00 1  IE  CO 
-1.0008E  or 
-1.0005E  <^0 
-1.0C03E  OC 
-1.0001E  OC 
-9.9987E-01 
-9.9973E-C1 
-9.9963E-C1 
-9.9956E-01 
-9.9954E-C1 
-9.9954E-C1 
-9.9957E-01 
-9.9961E-01 
-9.9967E-01 
-9.9974E-01 
-9.9981E-01 
-9.9987E-01 
-9.9993E-01 


Table  2.1.  (Continued) 
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K=50 


K=55 


K=A5 


K=36 


TIME 

CO 

l.OOOOE-Cl 
2.000OE-O1 
3.C000E-01 
4.0000E-01 
5.0000E-01 
6.0000E-01 
7.0000E-01 
P.OOOOE-Ol 
9.0CCCE-01 
1.0O00E  CO 
1.1COOE  00 
1.2000E  00 
1.3000E  00 
1.4000E  00 
1.5000E  00 
1.60C0E  00 
1.700CE  CO 
1.800CE  00 
1.9C00E  00 
2.0000E  00 
2.1000E  00 
2.2000E  00 
2.3000E  0" 
2.40C0E  00 
2.50C0E  00 
2.6C00E  00 
2.7000E  00 
2.8000E  00 
2.9000E  00 
3.0000E  00 
3.1000E  00 
3.2000E  00 
3.3000E  00 
3.40C0E  00 
3.5000E  00 
3.60G0E  CO 
3.7CC0E  00 
3.8C00E  ^0 
3.9000E  OC 
4.0000E  00 
4.1000E  00 
4.2000E  00 
4.3C00E  00 
4.4000E  00 
4.5000E  00 
4.6000E  00 
4.7000E  00 
4.80C0E  00 
4.9000E  00 
5.0000E  00 
5.1CCCE  00 
5.2000E  00 
5.300CE  00 
5.4000E  00 


OUT 
CO 

6.5182E-03 

4.1385E-02 

1.1165E-31 

2.1292E-01 

3.3623E-01 

4.7147E-01 

6.09UE-01 

7.4105E-01 

8.6100E-01 

9.6456E-01 

.0491E  30 

.1137E  00 

.1586E  00 

.1853E  00 

.1959E  00 

.1931E  00 

.1797E  00 

.1587E  00 

.1328E  00 

.1046E  00 

.0760E  00 

.0489E  00 

.0245E  00 

.0036E  00 

9.8672E-01 

9.7404E-01 

9.6544E-31 

9.6056E-01 

9.5895E-01 

9.6002E-01 

9.6318E-01 

9.6782E-01 

9.7336E-01 

9.7932E-01 

9.8525E-01 

9.9083E-01 

9.9581E-01 

l.OOOOE  00 

1.0034E  00 

1.0058E  00 

1.0C75E  00 

1.0084E  00 

1.0C86E  00 

1.0083E  00 

1.0075E  00 

1.0065E  00 

1.0053E  00 

1.0041E  00 

1.0028E  00 

1.0017E  00 

1.0007E  00 

9.9984E-01 

9.9918E-01 

9.9869E-01 


OUT 
0.0 

7.1698E-33 
4.5510E-02 
1.2269E-01 
2.3369E-31 
3.6830E-01 
5.1503E-01 
6.63C1E-01 
8.0304E-01 
9.2801E-01 
1.0331E  00 
1.1156E  00 
1.1747E  30 
1.2114E  00 
1.2276E  00 
1.2264E  00 
1.2112E  30 
1.1857E  00 
1.1536E  00 
1.1180E  00 
1.0820E  OC 
1.0478E  00 
1.0173E  00 
9.9169E-01 
9.7157E-01 
9.5715E-01 
9.4823E-31 
9.4430E-31 
9.4463E-01 
9.4837E-01 
9.5462E-01 
9.6250E-01 
9.7120E-01 
9.8C01E-01 
9.8836E-01 
9.9581E-01 
1.0021E  30 
1.0070E  00 
1.0105E  30 
1.0127E  00 
1.0136E  30 
1.0135E  00 
1.0126E  30 
1.0111E  00 
1.0C92E  00 
1.007CE  00 
1.0C49E  30 
1.0028E  00 
l.OOiOE  30 
9.9948E-01 
9.9828E-01 
9.9742E-31 
9.9689E-01 
9.9666E-01 
9.9669E-01 


OUT 

0.0 

5.8666E- 

03 

3.7258E- 

02 

1.0058E- 

01 

1.9206E- 

01 

3.0387E- 

■01 

4.2726E- 

■01 

5.5396E- 

■01 

6.7694E- 

•01 

7.9071E- 

•01 

8.9133E- 

01 

9.7634E- 

•01 

1.0446E 

00 

1.0961E 

00 

1.1317E 

00 

1.1529E 

30 

1.1617E 

00 

1.1601E 

30 

1.1505E 

30 

1.1350E 

30 

1.1158E 

00 

1.0946E 

00 

1.0729E 

30 

1.0520E 

30 

1.0328E 

30 

1.0159E 

30 

1.0018E 

00 

9.9052E- 

-01 

9.8215E- 

-01 

9.7648E- 

-01 

9.7324E- 

-01 

9.72C8E- 

-01 

9.7262E- 

-01 

9.7449E- 

-01 

9.7729E- 

-01 

9.8063E- 

-01 

9.8437E- 

-01 

9.8810E- 

-01 

9.9166E- 

-01 

9.9490E- 

-01 

9.9773E 

-01 

1.0001E 

00 

1.0C19E 

00 

1.0033E 

00 

L.0042E 

00 

1.0047E 

30 

1.0048E 

00 

1.0047E 

00 

1.0043E 

00 

1.0038E 

00 

1.0032E 

00 

1.0026E 

00 

1.0019E 

00 

1.0013E 

00 

1.0008E 

00 

OUT 

0.0 

4.6935E- 

•03 

2.9823E- 

02 

8.0610E- 

•02 

1.5427E- 

•01 

2.4494E- 

■01 

3.46C7E- 

•01 

4.5154E- 

•01 

5.5613E- 

•01 

6.5572E- 

■01 

7.4727E- 

•01 

8.2873E- 

•01 

8.9895E- 

•01 

9.5749E- 

■01 

1.0045E 

00 

1.0407E 

00 

1.0668E 

OC 

1.0842E 

00 

1.0941E 

00 

1.0977E 

00 

1.0963E 

00 

1.0911E 

00 

1.0832E 

00 

1.0736E 

00 

1.0630E 

30 

1.0522E 

00 

1.0416E 

00 

1.0316E 

00 

1.0226E 

3C 

1.0146E 

00 

1.0C78E 

30 

1.0023E 

00 

9.9790E- 

-01 

9.9459E- 

-01 

9.9227E- 

-01 

9.9081E- 

-01 

9.9007E- 

-01 

9.8992E- 

-01 

9.9023E- 

-01 

9.9089E- 

-01 

9.9179E- 

-01 

9.9283E- 

-01 

9.9394E- 

-01 

9.9506E- 

-01 

9.9613E 

-01 

9.9712E- 

-01 

9.9801E 

-01 

9.9878E- 

-01 

9.9943E 

-01 

9.9996E 

-01 

1.0004E 

00 

1.0007E 

30 

1.0C09E 

30 

1.0010E 

00 

1.0010E 

00 

Table  2.2.  Computer  output  data  for  K  =  55,  50,  45,  and  36 
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5.5000E  00 
5.6000E  00 
5.7000E  00 
5.8000E  00 
5.9000E  00 
6.0COOE  00 
6.1000E  00 
6.2000E  00 
6.3000E  00 
6.4000E  00 
6.5000E  00 
6.6000E  00 
6.700CE  00 
6.8000E  00 
6.9O00E  00 
7.0000E  00 
7.1000E  00 
7.2000E  00 
7.300CE  00 
7.4000E  00 
7.5000E  00 
7.60C0E  00 
7.7000E  00 
7.8000E  00 
7.9000E  00 
8.0000E  00 
8.1000E  00 
8.2000E  00 
8.3000E  00 
8.4000E  00 
8.5000E  00 
8.6000E  00 
8.7000E  00 
8.8C00E  00 
8.9C00E  00 
9.0000E  00 
9.1000E  00 
9.20C0E  00 
9.3000E  00 
9.4000E  00 
9.5000E  00 
9.6000E  00 
9.7000E  00 
9.8000E  00 
9.9000E  00 
l.OOOOE  01 


K=50 

OUT 
9.9838E-01 
9.9823E-01 
9.9821E-01 
9.9829E-01 
9.9846E-01 
9.9868E-01 
9.9893E-01 
9.9920E-01 
9.9945E-01 
9.9969E-01 
9.9989E-01 
1.0001E  00 
1.0002E  00 
1.0003E  00 
1.0003E  00 
1.0004E  00 
1.0004E  00 
1.0003E  00 
1.0003E  00 
1.0003E  00 
1.0002E  00 
1.0002E  00 
1.0001E  00 
1.000 IE  00 
l.OOOOE  00 
9.9998E-01 
9.9995E-01 
9.9994E-01 
9.9992E-01 
9.9992E-01 
9.9992E-01 
9.9993E-01 
9.9993E-01 
9.9994E-01 
9.9995E-01 
9.9997E-01 
9.9998E-01 
9.9999E-01 
9.9999E-01 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 


K=55 

9.9691E-01 

9.9729E-01 

9.9776E-01 

9.9828E-01 

9.9881E-31 

9.9931E-01 

9.9975E-01 

1.000 IE  00 

1.0004E 

1.0006E 

1.0008E 

1.0008E 

1.0008E 

1.0008E 

1.Q007E 

1.0005E 


00 
00 
00 
10 
00 
00 
00 
00 
00 
30 
00 
00 


1.0004E 

1.0003E 

1.0002E 

1.0001E 

9.9997E-01 

9.9989E-01 

9.9984E-01 

9.9981E-01 

9.9980E-01 

9.9980E-01 

9.9981E-01 

9.9984E-01 

9.9986E-01 

9.9989E-01 

9.9993E-01 

9.9996E-01 

9.9998E-01 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 

l.OOOOE  00 


K=45 

OUT 
1.0003E  00 
9.9992E-01 
9.9962E-01 
9.9940E-01 
9.9926E-01 
9.9919E-01 
9.9917E-01 
9.9920E-01 
9.9927E-01 
9.9936E-01 
9.9946E-01 
9.9957E-31 
9.9968E-01 
9.9978E-01 
9.9988E-01 
9.9996E-01 
l.OOOOE  00 
1.0001E  00 
1.0001E  00 
1.0C01E  00 
1.0001E  00 
1.0001E  00 
1.0001E  00 
1.0001E  00 
1.0001E  00 
1.0001E  00 
1.0C01E  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
9.9999E-01 
9.9998E-01 
9.9998E-01 
9.9997E-01 
9.9997E-01 
9.9997E-01 
9.9997E-01 
9.9998E-01 
9.9998E-01 
9.9998E-01 
9.9999E-01 
9.9999E-01 
9.9999E-01 
9.9999E-01 
l.OOOOE  00 


^=36 

OUT 
1.0010E  00 
1.0010E  00 
1.0009E  00 
1.0008E  00 
1.0007E  00 
1.0006E  00 
1.0005E  00 
1.0004E  00 
1.0003E  00 
1.0002E  00 
1.0001E  00 
l.OOOOE  00 
9.9998E-01 
9.9995E-01 
9.9992E-01 
9.9990E-01 
9.9989E-01 
9.9989E-01 
9.9989E-01 
9.9990E-01 
9.9991E-01 
9.9992E-01 
9.9993E-01 
9.9994E-01 
9.9995E-01 
9.9996E-01 
9.9997E-01 
9.9998E-01 
9.9999E-01 
9.9999E-01 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OOOOE  00 
l.OCOOE  00 
l.OOOOE  00 
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If  the  design  criteria  is  to  restrict  the  peak  overshoot  to 
+0.05  then  (2.9)  gives 

aK/K  '  HI?  '  14% 

Figure  2.5  shows  the  results  of  the  extreme  values  of  K  =  43 
and  K.  =  57  which  give  the  tolerable  -  limits  of  overshoot. 

Ill,   MULTILOOP  CONTROL  SYSTEM  SENSITIVITY 

A.  INTRODUCTION 

The  Kokotovic  [6]  method  of  sensitivity  analysis  is 
expanded  to  investigate  gam  and  time-contant  variations  in 
multiloop  systems  using  a  procedure  similar  to  the  one 
developed  for  the  basic  feedback  control  system.   This  proce- 
dure is  developed  using-.-the  block  diagram  of  a  regulated 
power  generator  which  incl-udes  an  amplifier-exciter,  a  series 
RC  compensating  circuit  in  the  feedback  loop,  and  a  load  as 
shown  in  Fig.  3.1. 

B.  GAIN-CONSTANT  SENSITIVITY  COEFFICIENTS 
Terms  are  defined  as  follows: 

Ki.Gi(s)  =   Direct  transfer  function  for  the  amplifier- 

exciter, 

KaG.zCs)'  =  Transfer  function  for  the  generator. 

K3G3(s)  =  Transfer  function  for  the  load. 

FH(s)  Feedback  c         >r  transfer  function. 

The  basic  equations  are  as  follows: 

c  K  v,  (s)KzG  !  s)K  G  (s) 

R~'~'      1+K:G  is;  FH(S)  -K,Gt  (s)K2G2  ^sy  K3G;<  (s;  (3.1a) 
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E2 

IT    (S) 

= 

=  1 

R"    (S) 

= 

Bi 

I+K1G1 (s)PH(s)+KiGi (s)K2G2 (s)K3G3 (s) 

I+K1G1 (s) FH(s) 
I+K1G1 (s)FH(s)+KiGi (s)K2G2 (s)K3G3 (s) 

K1G1  (s)FH(s) 
I+K1G1 (s)FH(s)+KiG! (s)K2G2 (s)K3G3 (s) 


(3.1b) 


(3.1c) 


(3. Id) 


As  before,  capital  letters  denote  the  Laplace  transform 
of  the  function  described,  except  for  K  and  F  which  are  gain 
constants,  and  T  which  is  a  time  constant. 

To  obtain  the  system  sensitivity  function  apply  (2.1) 
to  (3.1a)  to  give 

c  RK1G1K2G2K3G3 

1       (l+K1GiFH+K1G1K2G2K3G3) 2 

(3.2) 

C 

1+KiG!FH+KiGiK2G2K3G3 

This  sensitivity  function  corresponds  to  the  error  signal, 

E2,  in  the  original  multiloop  system  but  with  the  controlled 

output,  C,  as  the  input  to  the  system.   This  can  be  seen  by 

comparing  (3.2)  to  (3.1b)  with  C  replacing  R.   As  in  the 

basic  control  system  of  Section  II  ,  the  sensitivity  model  for 

gain  parameters  can  be  realized  by  cascading  two  identical 

systems  and  using  the  error  signals  in  the  sensitivity  model 

as  the  desired  gain  sensitivity  coefficients  as  shown  in 

Fig.  3.2. 

The  load-gain  sensitivity  function  is  derived  by  applying 

(2.1)  to  (3.1a),  yielding 
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RK1G1K2G2K3G3 (1+KiGiFH) 

T 

(3.3) 


SK3  (I+K1G1FH+K1G1K2G2K3G3) 2 


C(l+KiGiFH) 


I+K1G1FH+K1G1K2G2K3G3 


This  sensitivity  function  corresponds  to  the  error  signal, 
Ej,  as  can  be  seen  by  comparing  (3.3)  to  (3.1c)  with  C 
replacing  R  as  the  input. 

For  the  feedback-gain  sensitivity  function  apply  (2.1) 
to  (3.1a),  yielding 

RK1G1K2G2K3G3 (-K1G1FH) 
S    = 

I+K1G1FH+K1G1K2G2K3G3) 2 

(3.4) 
C(-KiGxFH) 


1+KiG1FH+K1GiK2G2K3G3 


This  sensitivity  function  corresponds  to  the  negative  of  the 
feedback  signal,  Bi,  as  can  be  seen  by  comparing  (3.4)  to 
(3. Id)  with  C  replacing  R  as  before. 

By  applying  (2.1)  to  the  system  transfer  function, 
sensitivity  functions  have  been  derived  for  the  various  gain 
parameters  and  the  pickoff  points  for  the  sensitivity 
functions  can  be  selected  by  inspection  as  shown  in  Fig.  3.2, 
No  new  components  have  been  required  and  all  sensitivity 
functions  can  be  generated  simultaneously. 
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C.   TIME-CONSTANT  SENSITIVITY  COEFFICIENTS 

Since  time-constant  variations  can  also  affect  the 
response,  sensitivity  coefficients  for  LR  and  CR  network 
combinations  can  add  helpful  information  to  the  analysis 
of  a  modeled  system.   For  this  derivation,  let 

C        KG         KG(Ts-t-l)  ,_   , 

R      1+KGFH      l+Ts(KGF+l)  U-bj 

Ts 
where  H  =  = — r-r     and  T  is  the  time  constant  of  an  RC  feedback 
Ts+1 

circuit.   Applying  (2.1)  to  (3.5)  yields 


S 


-R(KG) 2  FTs 
[1+Ts (KGF+1) ] 2 


C(-KGFTs) 


(3.6) 


[1+Ts (KGF+1) ]   Ts+1 

and  the  feedback-gain  sensitivity  is 

e       C(-KGFTs)  ,,  7, 

&F      1+Ts (KGF+1)  K  ' 

Therefore,  by  combining  (3.5)  and  (3.6) 

ST   =   SF  (TiTT»  (3-8) 

This  indicates  that  the  RC  time-constant  sensitivity  can 
be  generated  by  using  the  feedback-gain  sensitivity  as  the 
input  to  a  block  with  elementary  transfer  function  as  shown 
in  Fig .  3.2. 

In  the  direct  path  there  are  time-constants  associated 
with  the  generator  and  load  inductance  and  resistance 
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combinations.   The  effect  of  variations  of  these  time- 
constants  will  add  information  to  the  analysis  of  the  system, 
For  this  derivation,  let 

2   =     KG     =       K  (  r>    q\ 

R      1+KGFH      1+Ts+KFH  u   ' 

where  G  =     , ,  and  T  is  the  time-constant  of  a  direct  path 

LR  component.   Applying  (2.1)  to  (3.9)  yields 

-RKTs 
T      (1+Ts+KFH) 2 

C(-Ts) 
(1+Ts+KFH) 

and  the  direct-gain  sensitivity  is 


(3.10) 


s    =    C(Ts+l) (3  11) 

K       (1+Ts+KFH)  J.J-J-; 


Therefore,  by  combining  (3.10)  and  (3.11) 

S   (    ~Ts 
'T   ~    K  lTsTI 


**      =  SK  (b=It)  (3-12) 


As  with  the  feedback  RC  time-constant  sensitivity,  only  one 
elementary  block  has  to  be  added  in  series  with  the  direct- 
gain  sensitivity  model  to  complete  the  sensitivity  model  as 
shown  in  Fig.  3.2. 

D.   APPLICATION 

One  of  the  basic  problems  of  preliminary  design  analysis 
is  to  estimate  the  influence  of  each  design  parameter  on  the 
system  response  and  decide  how  improvement  of  the  system 
can  be  assured.   Sensitivity  coefficients  can  be  helpful  in 
solving  this  problem.   To  illustrate  this  procedure,  a  basic 
set  of  design  parameters  for  a  regulated  generator  is  taken  to 
be  amplifier  gam,  exciter  time-constant,  load-gain, 
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functions. 
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feedback-gain,  and  a  compensator  time-constant.   Figure  3.3 

is  a  computer  graph  of  the  time  response  to  a  unit-step 

input  and  the  five  design-parameter  sensitivity  coefficients. 

At  a  glance,  it  is  apparent  from  Fig.  3.3  that  Ki  and 

Ti  have  relatively  little  influence  on  the  response.   Table 

3.1  shows  the  maximum  value  of  S„   is  only  0.05  which  is  less 

than  10%  the  value  of  ST.  .   Figure  3.4  shows  Sv   and  S,_, 

K3  is-i        1 1 

where  the  scale  has  been  increased  by  10  to  show  the  general 
shape  of  the  curves.   Figure  3.5  shows  the  time  response  for 
Kx  =  40  and  Kj  =  60.   Even  for  this  50%  increase,  there  is 
little  change  in  the  performance.   Since  the  response  is 
very  "insensitive"  to  K1    and  Tx,    little  improvement  can  be 
achieved  by  adjusting  these  parameters  and  their  tolerances 
can  be  large  in  relation  to  the  other  parameters.   It  should 
be  noted  that  the  system  is  practically  insensitive  to  small 
variations  of  resistance  in  the  exciter  because  the  time- 
constant  and  gain  both  either  increase  or  decrease  by  the 
same  amount.   The  sensitivities  have  the  same  relative 
amplitude  but  differ  in  sign  and  therefore  the  effect  on  the 
time  response  is  canceled. 

In  Fig.  3.6  the  sensitivity  coefficients  for  the  load- 
gain  and  feedback-gain  constants  are  plotted.   Table  3.2 
shows  that  their  magnitudes  are  almost  equal  but  are  180 
degrees  out  of  phase.   Figure  3.7  is  a  graph  of  the  time 
response  with  their  original  values  and  the  response  for  a 
10%  increase  in  both  parameters.   Although  the  response  is 
sensitive  to  such  changes,  the  plots  are  almost  superimposed 
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due  to  the  one  being  the  negative  of  the  other.   This  could 
be  favorable  equalization  if  aging  or  temperature  cause 
these  gain  parameters  to  vary  by  the  same  percentage. 

Figure  3.3  shows  that  the  feedback-gain  and  the  compen- 
sator RC  time-constant  have  almost  the  same  sensitivity  but 
the  time-constant  sensitivity  function  is  delayed  by  one 
second. 

Since  fractional  parameter  variations  are  considered, 
the  permissible  percentage  of  a  parameter  change  can  be 
determined  from  the  allowable  system  response  variation. 

It  has  been  shown  that  any  general  form  of  a  control 
system  that  can  be  transformed  into  a  system  of  primary 
loops  can  be  analyzed  in  the  same  fashion  as  developed  in 
this  section.   The  sensitivity  to  variations  in  direct- 
and  feedback-gain  constants  are  available  as  error  or  feed- 
back signals  in  the  sensitivity  model  and  time-constant 
sensitivities  require  addition  of  elementary  transfer 
functions  as  shown  in  Fig.  3.2. 
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TIM 

0.0 

1.0000E-01 
2.0000E-01 
3.0000E-01 
4.0000E-01 
5.0000E-01 
6.0000E-01 
7.0000E-01 
8,0000E-0i 
9.0000E-01 
1.0000C  00 
1.1000E  00 
1.2000E  00 
1.3000E  00 
1.4000E  00 
1.5000E  00 
1.6000E  00 
1.7000E  00 
1.8000E  00 
1.9000E  00 
2.0000E  00 
2.1000E  00 
2.2000E  00 
2.3000E  00 
2.4000E  00 
2.5000E  00 
2.6000E  00 
2.7000E  00 
2.8000E  00 
2.9000E  00 
3.0000E  00 
3.1000E  00 
3.2000E  00 
3.3000E  00 
3.4000E  00 
3.5000E  00 
3.6000E  00 
3.7000E  00 
3.3000E  00 
3.9000E  00 
4.0000E  00 
4.1000E  00 
4.2000E  00 
4.3000E  00 
4.4000E  00 
4.5000E  00 
4.6000E  00 
4.7000E  00 
4.8000E  00 
4.9000E  00 
5.0000E  00 
5.1000E  00 
5.2000E  00 
5.3000E  00 
5.4000E    00 


Ski 


sTi 


0*0 

5.2083E-04 

2.9189E-03 

7.4409E-03 

1.4014E-02 

2.2527E-02 

3.2866E-02 

4.4918E-02 

5. 85706-02 

7.3713E-02 

9.0238E-02 

1.0804E-01 

1.2701E-01 

1.4705E-01 

1.6807E-01 

1.8995E-01 

2.1262E-01 

2.3597E-01 

2.5993E-01 

2.8440E-01 

3.0930E-01 

3.3455E-0X 

3.6007E-01 

3.8580E-01 

4.1166E-01 

4.3759E-01 

4.6351E-01 

4.8937E-01 

5.1512E-01 

5.4069E-01 

5.6603E-01 

5.9111E-0I 

6.1586E-01 

6.4025E-01 

6.6424E-01 

6.8779E-01 

7.1086E-01 

7.3344E-01 

7.5548E-01 

7.7696E-01 

7.9786E-01 

8.1816E-01 

8.3784E-0I 

8.5688E-01 

8.7528E-01 

8.9301E-01 

9.1006E-01 

9.2644E-01 

9.4214E-01 

9.5715E-01 

9.7147E-01 

9.8510E-01 

9.9804E-01 

1.0103E    00 

1.0219E    00 


0*0 

3.2895E-04 

1.2793E-03 

2.4687E-03 

3.7582E-03 

5.1151E-03 

6.5283E-03 

7.9903E-03 

9.4944E-03 

1.1034E-02 

1.2603E-O2 

1.4194E-02 

1.5801E-02 

1.7420E-02 

1.9043E-02 

2.0664E-02 

2.2280E-02 

2.38831-02 

2. 54 70 E -02 

2.7035E-02 

2.8574E-02 

3.0082E-02 

3.1556E-02 

3.2991E-02 

3.4385E-02 

3.5733E-02 

3.7032E-02 

3.8281E-02 

3.9476E-02 

4.0615E-02 

4.1696E-02 

4.2717E-02 

4.3677E-02 

4.4574E-02 

4.5408E-02 

4.6178E-02 

4.6881E-02 

4.7520E-02 

4.8093E-02 

4.8600E-02 

4.9042E-02 

4.9419E-02 

4.9731E-02 

4.9980E-02 

5.0168E-02 

5.0293E-02 

5.0357E-02 

5.0362E-02 

5.0309E-02 

5.0203E-02 

5.0039E-02 

4.9825E-02 

4.9561E-02 

4.9247E-02 

4.8887E-02 


0.0 

-3.2650E-04 
-1.  25781-03 
•2.4017E-03 
-3.6161E-03 
-4. 8670E-03 
•6.1426E-03 
•7.4349E-03 
■8.7368E-03 
-1.0042E-02 
•1.1343E-02 
•1.2635E-02 
-1.3910E-02 
■1.5165E-02 
■1.6393E-02 
-1.7590E-02 
•1.8751E-02 
•1.9872E-02 
■2.0948E-02 
■2. 1977E-02 
-2.2954E-02 
■2.3877E-02 
•2.4742E-02 
-2.5549E-02 
■2.6294E-02 
•2.6976E-02 
•2.7594E-02 
•2.8146E-02 
•2.8631E-02 
•2.9048E-02 
■2.9399E-02 
■2.9681E-02 
•2. 9896E-02 
•3.0044E-02 
■3.0126E-02 

3.0142E-02 
•3.0092E-02 
•2.9980E-02 
■2.9805E-02 
■2.9570E-02 

2.9276E-02 
■2.8925E-02 
•2.  8520E-02 

2.8062E-02 
■2.7553E-02 
■2.6997E-02 
•2.6393E-02 
■2.5746E-02 

2. 5059E-02 
■2.4335E-02 
•2.3572E-02 

2.2779E-02 
•2.1955E-02 

2. 1103E-02 
•2.0226E-02 


Table   3.1a.   Computer  output  data  for  C,   S^-,,   and    S^ 
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TIME' 

5.50006  00 
5.6000E  00 
5.7000E  00 
5.8000E  00 
5.9000E  00 
6.0000E  00 
6.1000E  00 
6.2000E  00 
6.3000E  00 
6.4000E  00 
6.5000E  00 
6.6000E  00 
6.7000E  00 
6.8000E  00 
6.9000E  00 
7.00006  00 
7.1000E  00 
7.2000E  00 
7.30006  00 
7.4000E  00 
7.5000E  00 
7.6000E  00 
7.70006  00 
7.8000E  00 
7.9000E  00 
8.0000E  00 
8.1000E  00 
8.2000E  00 
8.3000E  00 
8.4000E  00 
8.5000E  00 
8.6000E  00 
8.7000E  00 
8.8000E  00 
8.9000E  00 
9.0000E  00 
9.1000E  00 
9.2000E  00 
9.3000E  00 
9.4000E  00 
9.5000E  00 
9.60006  00 
9.7000E  00 
9.80006  00 
9.9000E  00 
1.0000E  01 
1.0100E  01 
1.02006  01 
1.030QE  01 
1.0400E 
1.05006 
►00 
00 
0800 
1.0900E  01 


1.061 
1.071 
1.081 


1.03286 
.04306 
.05266 
•0615E 
.06986 
.07756 
.08456 
•0910E 
.09696 
•1021E 
.10696 
.1111E 
.11486 
.11796 
.12066 
.1228E 
.12466 
.12596 
.12686 
.12736 
•1274E 
.12716 
•12666 
•1256E 
.12446 
.1229E 
•1211E 
.1191E 
•1168E 

..1143E 
•1116E 
.10876 
.10566 
.10246 
.09916 
.09566 
.09206 
.08836 
.0845E 
.08066 
.07676 
.0728E 
.06886 

..0647E 
.06076 
.05676 
.05276 
.0486E 

..0446E 
.0407E 
•0368E 
•0329E 
.02916 
.0253E 

1.0216E 


00 
00 
00 
00 

88 

00 
00 
00 
00 
00 
00 
00 
00 
00 

88 

00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 


*1 

4.8484E-02  - 


4.8035E-02 
4.7550E-02 
4 ,70241-02 
4.64676-02 
4.58756-02 
4.5251E-02 
4.4601E-02 
4.3923E-02 
4.3223E-02 
4.2501E-02 
4.1759E-02 
4.0999E-02 
4.0227E-02 
3.94436-02 
3.86496-02 
3. 78466-02 
3.70406-02 
3.62256-02 
3.54136-02 
3.4599E-02 
3.3787E-02 
3.2978E-02 
3.2178E-02 
3.13826-02 
3.05966-02 
2.98216-02 
2.90586-02 
2.83056-02 
2. 75696-02 
2.68 486-02 
2.6145E-02 
2.5457E-02 
2.4791E-02 
2.4144E-02 
2.3515E-02 
2.2910E-02 
2.2327E-02 
2.1766E-02 
2.12286-02 
2.07136-02 
2.02226-02 
1.9754E-02 
1.9313E-02 
1.88976-02 
1.8503E-02 
1.81346-02 
1.77926-02 
1.7472E-02 
1. 71796-02 
.69086-02 
•6664E-02 
•6441E-02 
.6244E-02 
1.6069E-02 


-{ 


aTi 

.93296-02 

•84086-02 

1.74756-02 

•1. 65236-02 

•1.55666-02 

1.45976-02 

•1.36206-02 

•1.26416-02 

■1.1660E-02 

•1.06816-02 

9.70376-03 

•8.73106-03 

•7.7655E-03 

•6.81066-03 

•5.86886-03 

•4.93921-03 

4.02446-03 

•3. 12856-03 

•2.24736-03 

1.38976-03 

•5.51036-04 

2.64216-04 

1.05646-03 

1.82136-03 

2.56236-03 

3.27456-03 

3.9590E-03 

4.6156E-03 

5.2448E-03 

5.8424E-03 

6.4102E-03 

6.9460E-03 

7.4535E-03 

7.9278E-03 

8.37066-03 

8.78556-03 

9.1657E-03 

9.5148E-03 

9. 83466-03 

1.0123E-02 

1.03826-02 

1.0610E-02 

1.08106-02 

1.09796-02 

1.11186-02 

1.12336-02 

1.13206-02 

1.13786-02 

1.14136-02 

1.14216-02 

1. 14076-02 

366E-02 

.3056-02 

.2216-02 

1.11166-02 


Table   3.1a.    (Continued) 
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TIME 

0.0 

1.0000 E-01 
2.0000E-01 
3.0000E-0 
4.0000E 
5.0000E 
6.0000E-01 
7.0000E-01 
8.0000E-01 
9.0000E-01 
l.OOOOE  00 
1.1000E  00 
1.2000E  00 
1.3000E  00 
1.4000E  00 
1.5000E  00 
1.6000E  00 
1.7000E  00 
1.8000E  00 
1.9000E  00 
2.0000E  00 
2.1000E  00 
2.2000E 
2.3000E 
2.4000E 
2.5000E 
2.6000E 
2.7000E  00 
2.8000E  00 
2.9000E 
3.0000E 
3.1000E 
3.2000E 
3.3000E 
3.4000E 
3.5000E 
3.6000E 
3.7000E 
3.8000E 
3.9000E 
4.0000E 
4.1000E  00 
4.2000E  00 
4.3000E  00 
4.4000E  00 
4.5000E  00 
4.6000E  00 
4.7000E  00 
4. 8000E  00 
4.9000E  00 
5.0000E  00 
5.1000E  00 
5.2000E  00 
5.3000E  00 
5.4000E  00 


0*0 

5.2081E-04 
2.9185E-03 
7.4408E-03 


-01  7.4408E-03  7 
-01  1.4014E-02  1 
-01      2.2527E-02      2 


00 
00 
00 
00 
00 


00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 


3.2866E-02 
4.4918E-02 
5.8570E-02 
7.3713E-02 
9.0238E-02 
1.0804E-01 
1.2701E-01 
1.4705E-01 
1.6807E-01 
1.8995E-01 
2.1262E-01 
2.3597E-01 
2.5993E-01 
2.8440E-01 
3.0930E-01 
3.3455E-01 
3.6007E-01 
3.8580E-01 
4.1166E-01 
4.3759E-01 
4.6351E-01 
4.8937E-01 
5.1512E-01 
5.4069E-01 
5.6603E-01 
5.9111E-01 
6.1586E-01 
6.4025E-01 
6.6424E-01 
6.8779E-01 
7.1086E-01 
7.3344E-01 
7.5548E-01 
7.7696E-01 
7.9786E-01 
8.1816E-01 
8.3784E-01 
8.5688E-01 
8.7528E-01 
8.9301E-01 
9.1006E-01 
9.2644E-01 
9.4214E-01 
9.5715E-01 
9.7147E-01 
9.8510E-01 
9.9804E-01 
1.0103E  00 
1.0219E    00 


sk5 

-0.0 
5. 2079 E -04 
2.9177E-03 
7.4349E-03 
.399QE-02 
.2462E-02 
3.2718E-02 
4.4627E-02 
5.8053E-02 
7.2863E-02 
8.8919E-02 
1.0609E-01 
1.2424E-01 
1.4323E-01 
1.6293E-01 
1.8323E-01 
2.0398E-01 
2.2507E-01 
2.4638E-01 
2.6780E-01 
2.8921E-01 
3.1050E-01 
3.3158E-01 
3.5235E-01 
3.7271E-01 
3.9257E-01 
4.1186E-01 
4.I051E-01 
4.4842E-01 
4.6555E-01 
4.8183E-01 
4.9721E-01 
5.1162E-01 
5.2505E-01 
5.3743E-01 
5.4874E-01 
5.5895E-01 
5.6803E-01 
5.7597E-01 
5.8275E-01 
5.8837E-01 
5.9281E-01 
5.9608E-01 
5.9819E-01 
5.9913E-01 
5.9892E-01 
5.9758E-01 
5.9513E-01 
5.9158E-01 
5.8696E-01 
5.8130E-01 
5.7463E-01 
5.6698E-01 
5.5839E-01 
5.4890E-01 


■1.9186E-04 
■1.6383E-03 
•4.9661E-03 


STH 

u.u 
-4.2725E-06 
-7.8408E-05 
-3.7232E-04 


:}:?§2if:8i  -i:S§Iii-8: 


2.6190E-0J 
-3.6636E-02 
-4. 8559E-02 
-6. 1829E-02 
-7.6317E-02 
-9. 1894E-02 
-1.0843E-01 
-1.2581E-01 
-1.4389E-01 
-1.6256E-01 
-1.8170E-01 
-2.0119E-01 
-2.2091E-01 
-2.4076E-01 
-2.6063E-01 
-2.8042E-01 
-3.0002E-01 
-3.1935E-01 
-3.3832E-01 
-3. 5684E-01 
-3.7483E-01 
-3. 9222E-01 
-4.0895E-01 
-4.2494E-01 
-4. 4014E-01 
-4.5449E-01 
-4.6795E-01 
-4. 8047E-01 
-4.9202E-01 
-5. 0256E-01 
-5.1206E-01 
-5.2051E-01 
-5.2788E-01 
-5.3415E-01 
-5. 3933E-01 
-5.4339E-01 
-5.4635E-01 
-5.4821E-01 
-5.4896E-01 
-5.4863E-01 
-5.4722E-01 
-5.4476E-01 
-5.4127E-01 
-5.3675E-01 
-5.3126E-0] 
-5.2480E- 
-5.1742E 
-5.0915E-01 
-5.0001E-01 


i-01 

;-oi 
i-oi 


-4.1039E-03 
-6.6987E-03 
-1.0113E-02 
-1.4404E-02 
-1.9608E-02 
-2.5750E-02 
-3.2838E-02 
-4.0866E-02 
-4.9819E-02 
-5.9671E-02 
-7.0385E-02 
-8. 1918E-02 
-9.4222E-02 
-1.0724E-01 
-1.2091E-01 
-1.3516E-01 
-1.4993E-01 
-1.6516E-01 
-1.8075E-01 
-1.9664E-01 
-2.1276E-01 
-2.2903E-01 
-2.4538E-01 
-2.6172E-01 
-2.7800E-01 
-2.9413E-01 
-3.1005E-01 
-3.2569E-01 
-3.4098E-01 
-3.5587E-01 
-3.7030E-01 
-3.8421E-01 
-3.9754E-01 
-4.1026E-01 
-4.2231E-01 
-4.3365E-01 
-4.4425E-01 
-4.5406E-01 
-4.6307E-01 
-4.7123E-01 
-4.7854E-01 
-4.8496E-01 
-4.9049E-01 
-4.9511E-01 
-4.9882E-01 
-5.0160E-01 
-5.0346E-01 
-5.0439E-01 
-5.0441E-01 


Table   3.1b.    Computer  output  data  for  0,    3^ ,    3f»   and   sTg 
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00 
00 
00 


TIME 

5.5000E  00 
5.6000E  00 
5.7000E  00 
5.8000E  00 
5.9000E  00 
6.0000E  00 
6.1000E  00 
6.2000E  00 
6.3000E  00 
6.4000E  00 
6.5000f 
6.60001 
6.70001 
6.8000E  00 
6.9000E  00 
7.0000E  00 
7.1000E  00 
7.2000E  00 
7.3000E  00 
7.4000E  00 
7.5000E  00 
7.6000E  00 
7.7000E  00 
7.  8000E  00 
7.9000E  00 
8.0000E  00 
8.1000E  00 
8.2000E  00 
8.3000E  00 
8.4000E  00 
8.5000E  00 
8.6000E  00 
8.7000E  00 
8.8000E  00 
8.9000E  00 
9.0000E  00 
9.1000E  00 
9.2000E  00 
9.3000E  00 
9.4000E  00 
9.5000E  00 
9.6000E  00 
9.7000E  00 
9.8000E  00 
9.9000E  00 
1.0000E  01 
1.0100E  01 
1.0200E  01 
1.0300E  01 
1.0400E  01 
1.0500E  01 
1.0600E  01 
1.0700E  01 
1.0800E  01 
1.0900E    01 


1.0328E  00 
1.0430E  00 
1.0526E  00 
1.0615E  00 
1.0698E  00 
1.0775E  00 
1.0845E  00 
1.0910E  00 
1.0969E  00 
1.1021E  00 
1.1069E  00 
1.1111E  00 
1.1148E  00 
1.1179E  00 
1.1206E  00 
1.1228E  00 
1.1246E  00 
1.1259E  00 
1.1268E  00 
1.1273E  00 
1.1274E  00 
1.1271E  00 
766E  00 
!56E  00 
►44E  00 
1.1229E  00 

bum  88 

1.1168E  00 
1.1143E  00 
1.U16E  00 
1.1087E  00 
1.1056E  00 
1.1024E  00 
1.0991E  00 
1.0956E  00 
1.0920E  00 
1.0883E  00 
1.0845E  00 
1.0806E  00 
1.0767E  00 
1.0728E  00 
1.0688E  00 
1.0647E  00 
1.0607E  00 
1.0567E  00 
1.0527E  00 
1.0486E  00 
1.0446E  00 
1.0407E  00 
1.0368E  00 
1.0329E  00 
1.0291E  00 
1.0253E  00 
1.0216E    00 


sk3 


5.3854E-01 

5.2736E-01 

5.1539E-01 

5.0269E-01 

4.8929E-01 

4.7525E-01 

4.6060E-01 

4.4540E-01 

4.2969E-01 

4.1352E-01 

3.9694E-01 

3.7999E-01 

3.6273E-01 

3.4520E-01 

3.2744E-01 

3.0951E-01 

2.9144E-01 

2.7329E-01 

2.5509E-01 

2.3690E-01 

2.1874E-01 

2.0066E-01 

1.8270E-01 

1.6490E-01 

I.4728E-01 

1.2990E-01 

1.1277E-01 

9.5940E-02 

7.9426E-02 

6.3260E-02 

4.7470E-02 

3.2081E-02 

1.7114E-02 

2.5921E-03 

-1.1465E-02 

-2.5040E-02 

-3.8116E-02 

-5.0678E-02 

-6.2716E^02 

-7.4216E-02 

-8.5169E-02 

-9.5567E-02 

-1.0540E-01 

-1.1468E-01 

-1.2338E-01 

-1.3151E-01 

-1.3907E-01 

-1.4606E-01 

-1.5247E-01 

-1.5833E-01 

-1.6362E-01 

-1.6836E-01 

-1.7255E-01 

-1.7620E-01 

-I.7932E-01 


•4.9006E-01 
■4.7932E-01 
•4.6784E-01 
4.5566E-01 
•4.4282E-01 
■4.2937E-01 
•4. 1535E-01 
•4.0079E-01 
3. 8576E-01 
3.7029E-0 
■3.5443E-0 
•3.3823E-0  . 
•3.2173E-0 
•3.  0497E-0 
•2. 8800E-0  . 
•2.7086E-0  . 
•2. 5360E-01 
•2.3625E-01 
■2.1887E-01 
•2.0148E-01 
•1.8414E-01 
•1.6687E-01 


1.49 


-Mi 


3272 
1.1590 
•9.9302E-02 
•8.2952E-02 
•6.6882E-02 
•5.1121E-02 
3.5691E-02 
2.0622E-02 
•5.9357E-03 
8.3437E-03 
2.2199E-02 
3.5609E-02 
4. 8554E-02 
6.1026E-02 
7.3006E-02 
8.4481E-02 
9.5444E-02 
1.0588E-01 
1.1579E-01 
1.2516E-01 
1.3399E-01 
1.4227E-01 
1.5001E-01 
1.5720E-01 
1.6385E-01 
1.6995E-01 
1.7551E-01 
1.8053E-01 
1.8502E-01 
1.8899E-01 
1.9245E-01 
1.9539E-01 


sTh 

5.0352E-01 
5.0172E-01 
4.9904E-01 
4.9549E-01 
4.9108E-01 
4.8584E-01 
4.7980E-01 
4.7296E-01 
4.6537E-01 
4.5705E-01 
4.4803E-01 
4.3834E-01 
4.2802E-01 
4.1710E-01 
4.0561E-01 
3.9359E-01 
3.8107E-01 
3.6810E-01 
3.5471E-01 
3.4094E-01 
3.2683E-01 
3.1242E-01 
2.9774E-01 
2.8283E-01 
2.6773E-01 
2.5247E-01 
2.371QE-01 
2.2165E-0I 
2.0616E-01 
1.9066E-01 
1.7518E-01 
1.5976E-01 


■1.4442E-01 

-1.2921E-01 

-1.1415E-01 

■9.9269E-02 

•8.4595E-02 

-7.0153E-02 

•5.5971E-02 

•4.2070E-02 

•2.8474E-02 

■1.5205E-02 

■2.2815E-03 

1.0278E-02 

2.2456E-02 

3.4237E-02 

4.5607E-02 

5.6552E-02 

6.7062E-02 

7.7126E-02 

8.6736E-02 

9.5883E-02 

1.0456E-01 

1.1277E-01 

1.2050E-01 


Table   3- lb-    (Continued) 
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IV.   EXAMPLES  OF  DIFFERENTIAL  EQUATION  SENSITIVITY  ANALYSIS 

A.   INTRODUCTION 

A  general  method  of  determining  the  sensitivity  functions 
for  the  coefficients  of  a  second-order  differential  equation 
is  derived  for  a  mass,  spring,  and  damper  system.   The 
procedure  is  then  applied  to  the  study  of  a  sixth-order 
differential  equation,  the  solution  of  which  is  initially 
unstable.   The  coefficients  are  then  adjusted  to  yield  a 
stable  solution  using  the  results  of  sensitivity  analysis. 

B".   SECOND-ORDER  DIFFERENTIAL  EQUATION  SENSITIVITY  TO 
COEFFICIENT  VARIATIONS 

The  elementary  mass-spring-damper  system  shown  in  Fig. 

4.1  is  described  by  the  following  differential  equation 

MX  +  BX  +  KX   =   u(t)  (4.1) 

Taking  the  Laplace  transform  and  solving  for  X(s)  gives 

X(s)  =     ^(S) (4.2) 

Ms  +Bs+K 

For  convenience,  let  M  =  1.   Applying  (2.1)  to  (4.2)  yields 

SR   =    =f  U(S)  2  =   X(s)  {-f^ )       (4.3) 

(s  +Bs+K)  S  +Bs+K 

By  comparing  (4.3)  to  (4.2)  it  can  be  seen  that  the  sensi- 
tivity model  is  a  replica  of  the  original  system  with  the 
controlled  output  as  the  input  to  the  model.   The  sensitivity 
coefficient  is  measured  at  the  point  corresponding  to  the 
feedback  signal. 
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U(s) 


Fig. 4.1.   Block  diagram   of  mass-spring-damper  system 
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Fig. 4. 2.    Second-order  differential  equation    sensitivity 
model. 
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Similarly, 

S    =   ~Bs  U(s) =   X(s)  (— ZSS )  (4.4) 

(s2+Bs+K)2  S2+Bs+K 

Figure  4.2  is  a  block  diagram  of  the  complete  sensitivity 

model  which  shows  clearly  that  S-.  and  SR  are  comparable  to 

the  feedback  gain  sensitivities  as  derived  in  Section  II 

and  III. 

The  solution  of  (4.1)  for  a  unit-impulse  input  is 

i        _  B         

x(t)   =   e    2  ^   sin  /~    B2  t         (4.5) 

rz — B2  K  "  T~ 

K    4~ 

Figure  4.4  is  the  graph  of  x(t)  and  the  sensitivity  coeffi- 
cients S„  and  S_.  for  B  =  0.6  and  K  =  0.6.   The  plot  shows 
that  the  system  is  more  sensitive  to  changes  in  K  than  B 
as  verified  by  (4.5).   Table  4.1  is  the  computer  output  data 
for  X/  SR,  and  Sg. 

C.   SIXTH-ORDER  DIFFERENTIAL  EQUATION  SENSITIVITY  TO  COEFFI- 
CIENT VARIATIONS 

Consider  the  differential  equation 

d6x     d5x     d"x     d3x     d2x      dx  /(4,6) 

a  +  b  +  c  +  d  +  e  +  f  —  +  g  x  =  u(t) 

dt6      dt5      dt"      dt3      dt2      dt 

where  nominally  a  =  1,  b  =  3.4,  c  =  11.4,  d  =  9.4,  e  =  7.4, 
f  =  0.575,  g  =  0.452.   The  solution  when  u(t)  is  a  unit-step 
can  be  determined  by  DSL  using  only  one  block.   For  the 
sensitivity  model  the  equation  must  be  integrated  in  six 
steps  to  isolate  each  sensitivity  coefficient.   The  procedure 
used  for  the  second-degree  equation  can  be  expanded  to  give 
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3.0 

2.0 

♦ 

X 

1.0 

0.0 

-1.0 

f«K 

IZ'jf 

SB 

,Jffv^_>^^^ 

JJS 

J20 
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-3.0 
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V 

/   SK 

-  ..  » 

TIME  (SEC) 

Fig. 4. 4.    Second-order  system  response   and    sensitivity 
functions. 
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Fig. A. 5.    Sixth-order  unstable    system  response   and 
sensitivity  functions. 

(a).   Response  and   sensitivity  functions  for 

d   and  f. 

(b).  Sensitivity  functions  for  b  and  o. 
(c).  Sensitivity  functions  for  e  and  g. 
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Fig. A. 5(a).  X,  Sd ,  and  Sf 
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Fig. 4. 6.   Sixth-order  stable    system  response   and 
sensitivity  functions. 

(a).   Response   and    sensitivity  functions  for 
d   and   f. 

(b).    Sensitivity  functions  for  b  and   c. 

(c).    Sensitivity  functions  for  e   and   g. 
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Fig. 4. 6 (a).   X,   Sd,   and   Sf 
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a  model  similar  to  Fig.  4.2.   The  expanded  sensitivity  model 
is  shown  in  Fig,  4.3.   The  solution  of  (4,6)  is  unstable  as 
can  be  seen  in  Fig.  4.5(a) ,   The  sensitivity  functions 
graphed  in  Fig.  4.5(b)  indicate  that  the  inner-loop  feed- 
back gain  parameters  b  and  c  influence  the  output  relatively 
little.   Going  from  the  coefficient  of  the  highest  order 
term  to  the  lowest,  the  influence  increases  but  it  is  ob- 
tained with  an  increasing  time  delay. 

The  problem  of  stabilizing  the  solution  may  be  solved 
using  sensitivity  analysis  by  noting  that  Sf  is  out-of -phase 
with  X.   By  increasing  the  coefficient  f,  the  response 
becomes  stable  as  shown  in  Fig,  4.6(a)  for  f  =  1.8.   Any 
further  adjustment  of  the  performance  can  be  achieved  by 
referring  to  the  sensitivity  functions  for  the  stabilized 
system  as  shown  in  Fig.  4.6. 
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TIMC         x  SENSE        SENSK 

0.0  CO         -0.0         -0.0 

1.0000^-01  4.8990E-03  -9. 6995E-05  -2.AA01E-06 

2.0000E-01  1.9185F-02  -7.5187F-C*  -3.8079E-05 

3.0000E-01  4.2229E-02  -2.A563E-03  -1.8789E-0A 

4.0000E-01  7.3386E-C2  -5.63C1*=-C3  -5.78A0E-0A 

5.00C0E-01  1.1200E-C1  -1.0622E-02  -1.37A6E-03 

6.0000E-01  1.57A2E-01  -1.7713E-C2  -2.7729E-03 

7.00C0E-0!  2.0898E-01  -2.7114F-02  -A.99A3E-03 

8.00C0F-01  2.6601F-C1  -3.8976C-C2  -8.2778E-03 

9.0000E-01  3.2788E-01  -5.3385E-02  -1.2874E-02 

l.OOCOE  00  3.9394E-C1  -7.0372E-C2  -1.9041E-02 

1.1000E  00  4.6356E-C1  -8.9912E-02  -2.7034E-02 

1.20C0E  00  5.3613E-C1  -1.1193F-C1  -3.7106E-02 

1.3000E  00  6.1106E-C1  -1.3632E-C1  -4.9500E-02 

1.A00CE  00  6.8776E-C1  -1.6290^-01  -6.4443E-02 

1.5000E  00  7.6570F-C1  -1.9150E-C1  -8.2148E-02 

1.6000^  00  8.AA34E-C1  -2.2187E-C1  -1.0280E-01 

1.70C0E  00  9.2319E-C1  -2.«5377E-01  -1.2657E-01 

l.nOOOE  00  1.0018E  00  -2.86Q2E-01  -1.5360E-01 

1.90C0E  OC  1.0797E  OC  -3.2102E-01  -1.8399E-01 

2.0000E  00  1.156AE  00  -3.5577P-C1  -2.1783E-01 

2.1000^  00  1.2317F  00  -3.9085E-C1  -2.5516E-01 

2.2000E  00  1.3051E  00  -4.259AE-01  -2.9600E-01 

2.3000E  00  1.3763E  00  -A.6072E-C1  -3.4034E-01 

2.4000E  00  1.A451E  CO  -4.9487E-01  -3.8812E-01 

2.5000E  00  1.5112E  00  -5.2807E-01  -4.3928E-01 

2.6000E  OC  1.57A3E  CO  -5.60C3F-01  -4.Q370E-0! 

2.70C0F  00  1.6343E  00  -5.9CA5E-01  -5.5123E-01 

2.8000E  00  1.6909F  OC  -6.19C5E-C1  -6.1173E-01 

2.9000E  00  1.7AA2E  OC  -6.4558^-Cl  -6.7498E-01 

3.0000F  00  1.7939F  CC  -6.6981E-01  -7.4077E-01 

3.1000E  00  1.839<3F  00  -6.91E1E-0I  -8.0885E-01 

3.2000E  00  1.8823E  00  -7.10A9F-01  -8.7898E-01 

3.3000E  00  1.92C9E  CC  -7.26«59»=-Cl  -9.5086F-01 

3.4000F  00  1.9558E  00  -7.3966E-01  -1.0242E  00 

3.5000E  OC  1.9870E  OC  -7.A959E-C1  -1.0987E  00 

3.6000E  00  2.0145E  CC  -7.5629E-C1  -1.1740E  00 

3.7000E  00  2.0384E  OC  -7.5968E-01  -1.2498F  00 

3.8000E  OC  2.0587E  CC  -7.5974E-01  -1.3258E  00 

3.9000F  00  2.0756E  CO  -7.56A6E-01  -1.4017E  00 

4.0000E  00  2.0891E  CC  -7.4985E-01  -1.4770E  00 

A.1000E  00  2.0994E  00  -7.3995F-C1  -1.5515E  00 

4.20CCE  00  2.1065F  CC  -7.2682E-C1  -1.6249E  00 

4.3000P  00  2.11C7E  CC  -7.1055E-C1  -1.6968E  00 

4.4000E  OC  2.1120E  CO  -6.9125E-01  -1.7669E  00 

4.5000E  00  2.11C7E  CC  -6.6906E-C1  -1.8349E  00 

4.60C0E  00  2.1068E  OC  -6.441CE-01  -1.9006E  00 

A.70C0E  00  2.1006E  00  -6.1656E-01  -1.9637E  00 

4.80C0F  00  2.0923E  00  -5.8661E-01  -2.0239E  00 

A.90C0E  00  2.082CE  CC  -5.5AA5F-01  -2.0809E  00 

«5.00CCC  00  2.06<?8E  00  -5.2027E-01  -2.13A7E  00 

5.1OC0E  00  2.0560E  00  -A.8A3CE-01  -2.18A9E  00 

5.20C0E  00  2.0AC8E  OC  -A.A676F-01  -2.2315E  00 

5.3OC0E  00  2.02A2E  00  -A.C788E-01  -2.27A2E  00 

5.AOC0E  00  2.0066E  00  -3.6789E-01  -2.3130E  00 


Table  4.1.  Computer  output  data  for  X,  Sb»  and  Sr 
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TIwp 

X 

SFNSR 

SFNSK 

5.5000*" 

OC 

1.988CE 

CC 

-3.27C3F-C1 

-2.347PE 

00 

5.60C0C 

00  : 

1.9686E 

cc 

-2„8!?55E-C1 

-2.3784= 

00 

5.7000E 

00   ! 

L.S485F 

cc 

-2.4366E-CI 

-2.A049E 

00 

5.80COE 

00   " 

1.928CF 

cc 

-2.0162F-C1 

-2.4271E 

00 

5.90C0F 

oc  : 

1.9C72E 

cc 

-1.5965F-01 

-2.4452F 

00 

6.000CF 

00   1 

:.R862F 

oc 

-1.1798F-C1 

-2.4591E 

00 

6.1000^ 

00 

L.8651F 

CO 

-7.6815«=-C2 

-2.4688E 

00 

6.200CF 

oc 

L.8441F 

OC 

-3.6382F-C2 

-2.4745E 

00 

6.300CE 

oc  : 

1.8233F 

00 

3.1268F-C^ 

-2.4761E 

00 

6.400CC 

oc 

I.R028E 

cc 

4.1C2CE-C2 

-2.4739E 

00 

6.5000E 

oc  : 

1.7827E 

cc 

7.8M6F-02 

-2.^67<5F 

00 

6.600CE 

oo  : 

L.7631E 

cc 

1.1425^-Cl 

-2.*582E 

00 

6.7000E 

oc 

1.7441F 

cc 

1.A827F-C1 

-2.4451E 

00 

6.8000E 

00 

L.7257E 

oc 

1.8C54F-C1 

-2.4296F 

00 

6.900CE 

oc 

..708CE 

cc 

2.1C92F-C1 

-2.40Q0E 

00 

7.00COE 

OC   1 

..6912«= 

oc 

2.3931F-C1 

-2.3865E 

00 

7.10C0E 

OC   ] 

..6752F 

CO 

2.65MF-C1 

-2.3612E 

00 

7.2000E 

oc  : 

..66CCE 

cc 

2.oc7«p>cl 

-2.3334E 

00 

7.300CE 

oc  1 

1.6458E 

cc 

3.11^3^-01 

-2.3033E 

00 

7.40CCE 

oc  1 

1.6326E 

cc 

3.3123E-C1 

-2.2712'= 

00 

7.5000E 

oo  : 

L.62C3F 

00 

3.4*51E-C1 

-2.2372c 

00 

7.60CCE 

oc 

L.6C89E 

cc 

3.63MF-C1 

-2.2016E 

00 

7.7000E 

oo  : 

..5986E 

CO 

3.76C1F-C1 

-2.1646E 

00 

7.8000E 

00 

..5e<;3F 

cc 

3.8f23r-01 

-2.1264c 

00 

7.90C0P 

oc  : 

L.581CF 

cc 

3.9^12F-C1 

-2.0874E 

00 

8.0000E 

00   1 

..5736F 

CO 

3.9972F-C1 

-2.0477E 

00 

8.10C0E 

oc 

L.5673F 

CC 

4.C3C6F-C1 

-2.0075E 

00 

8.20C0E 

00   1 

..C618E 

oc 

A.Ci2CF-Cl 

-1.9671E 

00 

8.30CCC 

OC   1 

..?5">4E 

CO 

4.0322P-C1 

-1.9268E 

00 

8.4000E 

00   , 

..C538E 

oc 

4.0C17E-01 

-1.P866E 

00 

8.5000E 

oc  ; 

L.551CE 

cc 

3.9516E-C1 

-1.8468E 

00 

8.60COE 

oo  : 

..5491E 

cc 

3.8828F-C1 

-1.8076F 

00 

8.70CCP 

00   i 

L.548GE 

cc 

3.7962F-.01 

-1.7692E 

00 

8.80CCE 

00   ! 

L.5477E 

cc 

5.6930E-C1 

-1.7317E 

00 

8.900CC 

00   1 

..5480F 

cc 

3.5743F-C1 

-1.6C54E 

00 

9.00C0F 

OC   1 

..5491E 

cc 

3.*4!3F-C1 

-1.6603E 

00 

Q.1000E 

oo  : 

L.55C7E 

cc 

3.2952F-C1 

-1.6266E 

00 

9.20CCE 

oc  : 

1.5529E 

00 

3.1373F-C1 

-1.59445 

00 

9.3000E 

00 

L.5557F 

cc 

2.9^89F-C1 

-1.S63QF 

00 

9.4000E 

00   1 

..5590E 

CO 

7.7913F-C1 

-1.5351E 

00 

9.50CCC 

oo  : 

L.  5627F" 

cc 

2.6C«r7F-Cl 

-1.5081F 

00 

9.6000E 

00   1 

..5667F 

oc 

2.413CE-C1 

-1.4830F 

00 

o.700CE 

oc 

L.5712F 

cc 

2.216CF-C1 

-1.4598E 

00 

9.80C0C 

oo  : 

L.5759E 

cc 

2.C144F-01 

-1.4387E 

00 

9.90CCE 

oc 

..5809E 

cc 

1.81CCC-C1 

-1.4105F 

00 

l.OOOOE 

01   1 

..5860F 

oc 

1.6C4CF-C1 

-1.^025= 

00 

1.0100E 

01 

..5014E 

CO 

1.3976F-C1 

-1.3875E 

00 

1.02C0F 

01 

..5969F 

CC 

1.192CE-C1 

-1.3745F 

00 

1.0300E 

01 

1.6024F 

CC 

9.8911F-02 

-1.3636E 

00 

1.0400E 

01 

..6081F 

CC 

7.8713E-02 

-1.35475 

00 

1.0500E 

01 

1.6137F 

cc 

5.8999F-02 

-1.3479E 

00 

1.0600E 

01   1 

L.6193E 

00 

3.9763E-C2 

-1.3429E 

00 

1.0700E 

01 

L.6249F 

oc 

2.1C90E-C2 

-1.3390E 

00 

1.0800* 

oi  : 

..6303E 

00 

3.0612F-C3 

-1.3387E 

00 

1.0900E 

01   1 

1.6357E 

cc 

-1.4252F-02 

-1.3392E 

00 

Table   4.1.    (Continued) 
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V.   NONLINEAR  SYSTEM  SENSITIVITY 

A.  INTRODUCTION 

While  it  is  true  that  no  physical  systems  are  exactly 
linear,  often  adequate  linear  approximations  can  be  used. 
In  these  cases,  the  Kokotovic  [6]  method  of  sensitivity 
analysis  as  developed  in  the  preceding  sections  of  this 
thesis  can  be  applied.   There  are  many  situations,  however, 
for  which  linear  representations  are  invalid  and  nonlinear 
characteristics  must  be  considered. 

Figure  5.1  is  a  block  diagram  of  a  nonlinear  system 
containing  a  linear  element  represented  by  a  transfer  function 

K  \    y    ,  where  D  =  d/dt  is  the  differential  operator,  and  a 

nonlinear  element,  N,  whose  output  is  denoted  by  f(e).   D  is 
used  instead  of  s  in  the  linear  transfer  function  because  the 
Laplace  transform  and  its  inverse  are  generally  not  suitable 
for  nonlinear  analysis. 

The  sensitivity  to  gain-constant  variations  is  developed 
for  the  general  nonlinear  system  and  then  applied  to  a 
specific  system  containing  a  saturation  element  as  the  non- 
linear component. 

B.  SENSITIVITY  COEFFICIENT  FOR  A  GAIN  FACTOR 

In  order  to  designate  the  dependence  of  the  variables  in 
Fig.  5.1  the  following  notation  is  used: 
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-    i\ 


Fig. 5.1.    Block,  diagram  of  nonlinear  system. 


»»    3K 


Fig. 5. 2.  Sensitivity  model  of  nonlinear  system 
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Let 


r  =  r(t) 

e  =  e(t,K)  =  r(t)  -  c(t,K) 

f(e)  =  output  of  nonlinear  element  N 

c  =  c(t,K)  =  r (t)  -  e(t,K) 

x  =  x(t,K) 


-.n         Hn_1  a 

P(D)   =   a  - —  +  a   , T   +  •••  +  ai^r  +  ao   (5.1a) 

n  ,.n    n-1  j.n-1  at 

dt         at 


Q(D)       =      bn^+Vi   g_L+    ...    +bl|_+b0  (5.1b) 


From  the  definition  of  gain-constant  sensitivity  coeffi- 
cient (2.1) 

bK   "   8K/K   "    *  8K  K°"C} 

Using  the  differential  equation  for  the  linear  element  of  the 
system  gives 

KP(D)  f(e(t,K))  =  Q(D)  c(t,K)  =  Q(D)  [r ( t) -e ( t ,K) ]      (5.4) 

Taking  the  logarithmic  derivative,  Z/d    in  K  =  K  d/%    K,  of 
(5.4)  and  substituting  (5.2)  yields 

KP(D)  f(e(t,K))  +  K2P(D)^fe(t'K)^ 
3e(t,K) 


3K 

(5.5) 
-KQ(D) 


Rearranging  and  simplifying  (5.5)  gives 

q        _^l£      KP(D)    ,      K2P(D)  3f  (e)  (5.6) 

K        K8K   "    Q(D)  Z[e)     +      Q(D)   9K 

Figure  5.2  is  a  block  diagram  of  the  proposed  sensitivity 

model  for  the  nonlinear  system.   It  will  be  shown  that  the 

describing  equation  of  this  model  is  identical  to  (5.6)  . 
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From  Fig.    5  .  2 

sk    "    Wic    "    °-b  <5-7' 

Substituting  (5.2)  in  (5.7)  yields 

b  =  c  +  K  ~  (5.8) 

Applying  the  transfer  function  for  the  linear  element  of 
the  model  gives 

-K2P(D)  !|i^!  ||  =  Q(D)  b  =  Q(D)c  +  KQ(D)  ||  (5.9) 

Substituting  (5.4)  in  (5.9)  for  Q(D)c  and  solving  for  S^ 
as  defined  by  (5.2)  gives 

SK   ~     "^K  "  K  qTdT   (S)    K  0W  3K  (5*10) 

Equations  (5,10)  and  (5.6)  are  the  same  and,  hence,  the  model 
generates  the  gain-constant  sensitivity  coefficient  of  inter- 
est.  The  sensitivity  model  block  represented  by  « — - — -   is  the 
instantaneous  slope  of  the  nonlinear  function  f(e).   In  many 
nonlinear  systems  a  piecewise-linear  approximation  can  be 
utilized.   The  sensitivity  model  will  be  the  same  as  that  of 
a  linear  system  except  that  the  gain  changes  as  a  function  of 
the  error  signal  in  the  original  nonlinear  system.   An  example 
for  a  piecewise-linear  saturable  element  will  make  this  more 
meaningful. 

C.   APPLICATION 

The  nonlinear  component  used  as  an  example  will  be  a 
limiter  with  a  saturation  characteristic  as  shown  in  Fig.  5.3. 
The  linear  plant  used  in  Section  II  is  used  for  comparison 
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of  the  time  responses  and  gain  sensitivity  coefficients. 

The  transfer  function  in  differential  operator  notation 

becomes 

KP  (D)      50 

Q(D)       D3  +  10D2  +  24D 

The  gain-constant  sensitivity  model  as  shown  in  Fig.  5.2 
requires  a  piecewise-linear  function  described  by  Fig.  5.4. 

It  can  be  seen  that  -s — is  the  slope  of  the  nonlinear 

function  f (e)  and  may  be  handled  by  logic  statements  in  the 
DSL  computer  program  as  given  in  Appendix  D. 

Figure  5.5  shows  the  time  response  and  the  sensitivity 
function  for  a  unit-step  input  with  a  limited  amplitude  of 
error  signal  set  at  0.1.   The  nonlinear  element  is  saturat- 
ed until  4.8  seconds.   During  this  period  of  saturation  the 
gain  representing  the  nonlinear  element  in  the  sensitivity 
model  is  zero,  and  therefore,  the  sensitivity  coefficient 
and  the  output  of  the  nonlinear  system  are  one  and  the  same 
as  shown  in  Fig.  5.5  for  the  first  4.8  seconds.   Figure  5.6 
shows  the  response  and  sensitivity  for  the  limited  amplitude 
of  pi  =  -1.0  and  pz    =  1.0.   This  lets  the  system  perform  in 
a  completely  linear  fashion  and  the  effect  of  the  nonlinear 
element  may  be  realized  by  comparison  with  Fig.  5.5.   The 
nonlinear  response  is  slowed  down,  the  maximum  sensitivity 
to  gain  is  greater  but  delayed  by  3.8  seconds;  also,  the 
peak  overshoot  is  greatly  reduced. 

To  verify  the  sensitivity  coefficients  for  gain  vari- 
ations, check  points  were  chosen  in  the  saturation  region, 
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f(e)    n 


f(e)=px,  e<px 
f(e)=p2,  e>P2 
f(e)=e,        Pi^e<p2 


?i8«5.3.   Characteristics   of  nonlinear  N. 


9f(e)    i 

i 

de 

1 

e 

i 

Pi 

P2 

3f(e) 
?e 

*f(e) 
de 

€XlgJ 

ae 


1 ,      P]£e  <  P2 
0 ,      e  <  vi 
0 ,      e  >  p2 


Fig. 5. 4.   Characteristics  of  sensitivity  model  nonlinearity. 
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1.5 


1.0 


0.5 


0.0 


^ 

/SK 

- 

?02 

)0<                      V 

106          / 

m    t(sec) 

-0.5 

Fi6«5.5.   Response  and   sensitivity  for  nonlinear  system 
with  pi=-0.1  and   p2=  0.1. 
1.5 


1.0 


0.5 


0.0 


-0.5 

Fig. 5.6.  Response  and  sensitivity  for  unlimited  linear 
system. 
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peak-positive,  and  peak-negative  portions  of  the  sensitivity 

curve  as  shown  in  Fig.  5.7  for  pi  =  -0.2  and  p2  =  0.2.   The 

approximate  change  in  the  response  due  to  a  fractional  change 

of  0.1  in  the  gain  parameter, 

AC  -  S^  AK/K 
is. 

was  calculated  and  compared  to  the  computer  output  for  K  =  55 

■ 

and  K  =  4  5  as  given  in  the  Table  5.2.   The  maximum  error 
found  was  0.5%  for  point  2  as  presented  in  Table  5.1. 

Figure  5.8  shows  the  response  for  14%  variations  in 
gain  as  was  required  for  the  linear  system  to  constrain  the 
peak  overshoot  to  within  +  0.05.   The  peak  overshoot  for 
the  nonlinear  system  was  constrained  to  +  0.014  for  the  same 
gain  variations.   This  can  be  predicted  from  Fig.  5.8  since 
at  the  time  of  peak  overshoot,  the  sensitivity  function  is 
much  nearer  zero  than  the  sensitivity  function  for  the  linear 
system  shown  in  Fig.  5.7. 
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1.5 


1.0 


0.5 


0.0 


-0.5 


^C 

12. 

^H 

J 1 

m 

zzi            \ 

3 

set 

3CB 

TIMS (SEC) 


Fig.   5 *7.   response   and    sensitivity  for  nonlinear  system 
with  p]_=-0.2   and   P2=   0.2. 


Pt. 

T 

G50 

sK 

C+AC 

G55 

C-AC 

C45 

1 

1.5 

.45215 

.45215 

.49737 

.49736 

.40694 

.40693 

2 

2.5 

.86763 

.84752 

.95238 

.94755 

.78288 

.78127 

3 

4.0 

1.0177 

-.21111 

.99659 

.99673 

1.0388 

1.0351 

Table   5.1.   Comparison  of  actual  outputs  and   outputs 
predicted   from   sensitivity  coefficients. 
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1.2 


1.0 


0.8 


0.6 


0.4 


0.2 


CO 


4  6 

TIME  ( S£C ) 


Fig. 5. 8.   Gain  tolerance   of  14$  for  overshoot  of  ±0.014. 


62 


K=50  K=55  K=45 


TIMF 

SFNSK 

mi 

cut 

OUT 

0,0 

CO 

0.0 

CC 

CC 

1.0000E- 

•01 

1.30ME-C3 

1.3CA1E- 

■C3 

1.A3A5F- 

-C3 

1.1737F- 

•03 

2.00C0E- 

•01 

8.302AE-C? 

8.3024F- 

•03 

9.1326F- 

•03 

7.A721E- 

•03 

3.0000E- 

•01 

2.2552F-C2 

2.2552F- 

•C2 

2.APC7P- 

-C2 

2.C297E- 

•C2 

A.OOCOF- 

•01 

4.3547E-02 

4.35A7F- 

-C2 

4.7902P- 

-02 

3.9193F- 

•C2 

F.OOOOE- 

•01 

7.0099E-02 

7.0C99E- 

•02 

7.71C9F- 

-02 

6.3089F- 

•02 

6.0000E- 

•01 

1.0094E-01 

1.0C94F- 

-CI 

1.11CAE- 

-CI 

9.08APE- 

■02 

7.0000F- 

•01 

1.3498E-C1 

1.3A98E- 

■CI 

1.A8A7F- 

-01 

1.21A8F- 

■01 

8.0000E- 

•01 

1.7132E-C1 

1.7132E- 

■01 

1.88A5F- 

-01 

1.5A19E- 

•01 

9.0000E- 

•01 

2.0930E-C1 

2.093CC- 

■CI 

2.3023F- 

-CI 

1.8837F- 

■01 

l.OOOOE 

00 

2.48A3F-C1 

2.A8A3F- 

•01 

2.7328F- 

•CI 

2.235CE- 

•CI 

1.1000F 

00 

2.8P37E-C1 

2.8837E- 

•CI 

3.1721E- 

•01 

2.5953c- 

■CI 

1.2000E 

oc 

3.2886E-C1 

3.2886F- 

-01 

3.6174E- 

-01 

2.9597E- 

•01 

1.3000E 

oc 

3.6972E-C1 

3.6972c- 

•CI 

A.067CF- 

-CI 

3.3275F- 

•C! 

l.AO00c 

00 

A.1085F-C1 

A.1085F- 

■01 

A.5193E- 

-01 

3.6976F- 

•CI 

1.5000F 

00 

4.5215E-01 

4.5215E- 

•CI 

A.9736F- 

-01 

A.C693P- 

•01 

1.60C0E 

00 

4.9357E-C1 

A.9357F- 

•CI 

5.A292F- 

-01 

A.AA21F- 

•CI 

1.7000E 

00 

5.3507E-C1 

5.35C7F- 

•CI 

5.8857F- 

-CI 

A.9156E- 

•C! 

1.8000E 

00 

5.7662E-01 

5.7662F- 

-01 

6.3A28F- 

-01 

C1896F- 

•01 

1.9000F 

00 

6.1821E-C1 

6.1821F- 

-CI 

6.80C3E- 

-CI 

5.,5639F- 

•01 

2.0000E 

00 

6.5983E-C1 

6.5983F- 

-CI 

7.2*81F- 

-CI 

5.S385P- 

•01 

2.1000E 

00 

7.0146E-C1 

7.C1A6F- 

-CI 

7.7161F- 

-CI 

6.3132F- 

•CI 

2.2000E 

00 

7.A311E-C1 

7.A311F- 

-CI 

8.17A2F- 

-CI 

6.688CE- 

•CI 

2.30COE 

CO 

7.8476F-C1 

7.8^76F- 

•CI 

8.6294E- 

•01 

7.062eF- 

•CI 

2.A000E 

oc 

8.2A9AE-01 

8.26^1F- 

-CI 

9.C691E- 

■CI 

7.A377F- 

-CI 

2.500CE 

CO 

8.4752E-CI 

8.6763E- 

-CI 

9.A755F- 

-01 

7.8127F- 

-01 

2.6000E 

CO 

P.3906E-01 

9.C72CC- 

-CI 

9.833CF- 

-CI 

8.1877F- 

-CI 

2.7000c 

oc 

7.9771«=-C1 

9.A366F- 

-01 

1.013C* 

CO 

8.56CCE- 

-01 

2.8O00E 

oc 

7.2738E-C! 

9.7577F- 

-CI 

1.0361E 

CC 

8.92C5E- 

-CI 

2.9000E 

oc 

6.3A57E-C1 

1.0C27F 

CO 

1.C523F 

CC 

9.2573F- 

•CI 

3.0000E 

00 

5.2669E-C1 

1.C239F 

CO 

1.C621F 

CO 

9.5599F- 

•CI 

3.100CE 

00 

A.1113E-01 

1.0395E 

CC 

1.C661F 

CO 

9.82C7F- 

■01 

3.2O00E 

00 

2.9A68E-01 

1.0A96F 

CO 

1.0652F 

CO 

1.0035F 

oc 

3.3000E 

oc 

1.8319F-01 

l.C«5*PF 

CC 

1.06C3F 

CC 

1.02C3F 

00 

3.A000F 

CO 

fl.13A4F-C2 

1.05C8P 

CC 

1.0526E 

CO 

1.C32AE 

CC 

3.500CF 

CO 

-7.3A81E-C3 

1.C533F 

CC 

1.0A31F 

CC 

1.0AC3F 

CC 

3.600CE 

00 

-8.0616F-C2 

1.0A83E 

CC 

1.032PP 

CC 

l.C^AAE 

CC 

3.7000F 

00 

-1.3732E-C1 

1.CA16F 

CO 

1.C22AF 

CC 

1.CA53E 

CC 

3.8000P 

00 

-1.7736F-C1 

1.C33PE 

00 

1.C126F 

CO 

1.CA36F 

CO 

3.9000E 

oc 

-2.01A7E-C1 

1.0256F 

CO 

l.OCAC^ 

00 

1.0AC1E 

CO 

A.OOCOE 

00 

-2.1111E-C1 

1.0177F 

CO 

9.9673F- 

-01 

1.C351E 

CC 

A.1000E 

oc 

-2.0823E-C1 

1.C1C3F 

oc 

9.9112F- 

-01 

1.C29AF 

CO 

4.2000E 

00 

-1.95C9F-C1 

1.CC3PF 

CO 

9.8717F- 

-01 

1.C23AP 

CO 

A.3000E 

00 

-1.7412F-01 

9.98A5F- 

-CI 

9.PA79E- 

■CI 

l.C17^»E 

CC 

A.AOOOF 

oc 

-1.4774E-C1 

9.9A25E- 

-CI 

9.8382F- 

-CI 

1.C118F 

00 

A.SOOOE 

00 

-1.1826E-C1 

9.912AF- 

-CI 

9.8AC6F- 

-CI 

1.0067E 

CC 

A.60CCE 

00 

-8.7735E-02 

9.893AF- 

-CI 

9.8526F- 

-01 

1.0C2AF 

CO 

4.7000E 

00 

-5.79*53^-02 

9.88A4F- 

•01 

9.871AF- 

-01 

9.9881F- 

-CI 

A.800CF 

00 

-3.03A5E-C2 

9.883<5F- 

•CI 

9.99A7E- 

-01 

9.96C7E- 

-CI 

4.9000F 

00 

-5.9862E-03 

9.89C1F- 

-01 

9.92CCF- 

-01 

9.9A12E- 

-01 

S.OOCOP 

00 

1.AA02E-02 

9.9014E- 

-01 

9.9A5AF- 

•CI 

9.9289E- 

-CI 

5.10C0E 

oc 

3.0A35E-C2 

9.9161F- 

-01 

9.9693E- 

•01 

9.9229E- 

-01 

5.2OC0E 

00 

4.2025E-C2 

9.9327E- 

-CI 

9.99C5F- 

-CI 

9.9222F- 

-01 

5.3000E 

00 

4.9333E-C2 

9.9A97F- 

-01 

1.00C8F 

CO 

9.9257F- 

-CI 

5.A0C0E 

00 

5.2723F-02 

9.9662F- 

-01 

1.0022F 

00 

9.9323F- 

-CI 

Table   5.2.   Computer  output  data  for  nonlinear  system 
outputs  and    sensitivity. 
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TIMp 

5.f00CE  00 

5.60CCE  00 

5.70C0E  00 

5.80CCE  00 

•S.90C0E  00 

6.0000E  00 

6.100CE  00 

6.200CE  00 

6.3000E  00 

•S.AOOOE  00 

6.5000E  00 

6.600CE  00 

6.7000E  00 

6.8000E  00 

6.9000E  00 

7.0000E  00 

7.1000E  00 

7.20C0E  OC 

7.3000E  00 

7.4000E  00 

7.5000E  00 

7.6O00E  00 

7.7000E  00 

7.8000F  00 

7.90C0E  00 

8.00CCE  00 

8.1OC0E  OC 

8.2OC0P  00 

8.30C0E  00 

8.A00CE  00 

8.5000F  00 

8.60C0E  00 

8.7OC0E  00 

8.80COE  00 

8.9000E  00 

9.0000E  00 

9.1000E  00 

9.2OC0E  OC 

9.3000E  OC 

Q.AOCOF  OC 

9.5000E  00 

9.60C0E  00 

9.7000E  00 

9.8000^  00 

9.90C0E  00 

l.OOOOE  01 

1.01 OOE  01 

1.0200E  01 

1.03C0E  01 

1.0400F  01 

1.0500E  01 

1.0600P  01 

1.0700E  01 

1.0800E  01 

1.0900E  01 


SENSK 

5.27C3F-C2 

4.9877E-C2 

4.4889E-C2 

3.8392E-C2 

3.1CC1F-C2 

2.3272E-C2 

1.5682E-C2 

8.61A1E-C3 

2.3564E-C3 

-2.8982E-C3 

-7.04A8F-03 

-1.0060E-02 

-1.1984E-C2 

-1.2908E-02 

-1.2962E-C2 

-1.2303E-C2 

-1.1C94E-C2 

-9.5025E-03 

-7.6836F-C3 

-5.7786E-03 

-3.9067E-03 

-2.1646E-C3 

-6.2460E-04 

6.6608E-C" 

1.6829E-C3 

2.4202E-C3 

2.8889F-C3 

3.112AE-C3 

3.1238E-03 

2.9615E-C3 

2.6658E-C3 

2.2785E-C3 

1.8387E-C3 

1.3797E-C3 

9.3007E-C* 

5.1326E-04 

1.4567E-CA 

-1.6022E-04 

-4.0054E-C* 

-•5.7411E-C4 

-6.8283E-04 

-7.3338F-C* 

-7.3433E-C4 

-6.9433E-G4 

-6.2269F-0* 

-5.29e9E-C4 

-4.2456F-C* 

-3.164AE-C4 

-2.1058E-C* 

-1.1283E-C* 

-2.7299E-G5 

4.3333E-05 

9.8586E-C5 

1.3816E-C4 

1.628AE-C4 


K=50 

OIT 
9.9P13F-01 
9.QQ44F-01 
1.00C5E  CC 
1.0013E  CC 
1.CC19F  CO 
1.0023E  CC 
1.0024E  CC 
1.0C2AF  00 
1.0C23F  CC 
1.002CE  00 
1.CC17E  OC 
1.0013E  CO 
1.0C1CF  OC 
1.00C6F  00 
1.00C3E  OC 
1.0001F  CC 
9.9985F-C1 
9.9969P-01 
9.99e5PF-Cl 
9.9951E-C1 
9.99A9E-01 
9.995CE-C1 
9.9953F-C1 
9.99*9F-01 
9.9966F-01 
9.9973F-C1 
9.9981F-C1 
9.9<3flf?F-ci 
9.9994F-C1 
9.9999E-01 
l.OOCCE  00 
1.0CC1F  CC 
1.0CC1F  00 
i.OOOlF  CC 
1.0CC1F  00 
1.CCC1F  CC 
1.0C01E  00 
I.OOOIE  CC 
1.0CC1F  CO 
1.0GC1E  CC 
l.OCOCE  00 
l.OOCCF  CC 
l.OOOOF  00 
l.OCCCF  00 
9.9999F-C1 
9.9998F-C1 
9.9998F-C1 
9.999eP-01 
9.9C<58F-01 
9.9998E-01 
9.999PE-01 
9.9998E-01 
9.9998E-C1 
9.9999E-C1 
9.9999F-C1 


K=55 

OIT 
1.0031F  CC 
L0037E  00 
1.C0ACE  CO 
1.0C39F  CC 
1.0036F  CO 
1.CC31F  00 
1.C026F  CO 
1.0C19F  00 
1.0013E  00 
1.0CC7F  OC 
1.00C2E  CO 
9.998CF-C1 
9.9946F-01 
9.9Q23E-01 
9.99C9E-C1 
9.99C3F-01 
9.9904F-C1 
9.9912E-01 
9.9923F-01 
9.9937^-Cl 
9.9952F-C1 
9.9C67E-01 
9.9982F-C1 
9.9994F-C1 
l.OOOCE  00 
l.OCCl*  CO 
1.0CC2P  CO 
1.0CC2E  CO 
1.00C2E  CO 
1.C0C2E  00 
1.0CC2E  CO 
1.CCC2F  00 
1.00C2F  CO 
1.00C1F  00 
I.OOOIE  00 
l.COCCE  CO 
l.COCCP  CO 
9.9999E-01 
9.9997E-01 
9.9Q95E-C1 
9.999AE-C1 
Q.9994F-01 
9.Q994F-01 
9.999A»=-C1 
9.9995F-C1 
9.9996E-C1 
9.9997E-C1 
9.9998E-01 
9.9999F-01 
9.99C9E-C1 
1.000C«=  00 
l.OOOCE  CO 
l.OOOCE  CO 
l.OOCCF  CO 
l.OOOCE  CO 


K=45 

OIT 
9.9411F-C1 
9.9511F-01 
9.9616P-C1 
9.9718E-C1 
9.9813E-C1 
9.9899E-01 
9.9971F-C1 
1.00C3E  CO 
1.00C7E  00 
1.0011E  CC 
1.0012E  00 
1.0013F  00 
1.0013F  00 
1.CC13P  CO 
1.CC11F  00 
1.001CF  CC 
1.C0C8F  00 
1.00C6F  CO 
1.00C5F  00 
1.00C3F  CC 
1.C0C1F  00 
l.OOCCE  00 
9.9993F-01 
9.9986E-C1 
9.9981E-01 
9.9978F-01 
9.9977F-C1 
9.9977E-C1 
9.9978E-01 
9.998CF-C1 
9.9983E-C1 
9.9986F-C1 
9.9989F-C1 
9.9992F-01 
9.9995E-01 
9.9997E-01 
9.9999F-C1 
l.COOCE  00 
l.OOOCE  CO 
l.OCCCF  CO 
1.00CCC  CO 
l.OCCCF  CC 
l.OOCCE  CO 
l.CCOCE  CO 
l.OOOCF  CO 
l.OCCCF  CO 
l.OCOCF  CC 
l.CCCCE  00 
l.OOOCF  CO 
l.CCCCF  CO 
l.OOOCF  CO 
l.OCOCE  00 
9.9999E-01 
9.9999F-01 
9.9999E-01 


Table  5.2.  (Continued) 
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VI.   CONCLUSIONS 

The  technique  of  generating  sensitivity  functions  for 
a  basic  linear  feedback  control  system  using  the  concept 
of  a  sensitivity  model  is  developed.   This  model -is  a  replica 
of  the  original  system  and  is  cascaded  with  it.   All  sensi- 
tivity functions  are  generated  simultaneously  with  the  time 
response  by  use  of  the  IBM/360  Digital  Simulation  Program. 

Sensitivity  of  a  differential  equation  with  respect  to 
its  coefficients  is  investigated.   Computer  graphs  show  the 
differences  in  the  influence  of  the  coefficients  on  the 
solution  of  the  differential  equation.   The  coefficients  of 
the  lower-order  terms  have  the  greatest  influence  but  the 
effect  is  obtained  with  the  greatest  time  delay. 

A  theoretical  sensitivity  model  is  developed  for  a  non- 
linear control  system.   Instead  of  a  replica  of  the  original 
control  system  as  required  for  the  linear  model,  the  non- 
linear element  is  replaced  by  a  component  with  the  instanta- 
neous gain  equal  the  slope  of  the  nonlinear  function.   An 
example  illustrates  the  procedure  for  a  nonlinearity  due  to 
a  saturation  element. 

Sensitivity  analysis  indicates  that  proper  interpretation 
of  the  sensitivity  function  results  in  accurate  estimates  for 
changes  in  the  system  response  due  to  parameter  variations. 
Tolerance  limits  of  parameter  variations  are  determined  and 
performance  improved  by  application  of  sensitivity  coefficients 
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APPENDIX  A 


DSL/360  DIGITAL  SIMULATION  PRHGRAM  FOR  BASIC  SYSTEM 

TITLE  BASIC  SEPVC  SENSITIVITY  CGECFICIENTS 

INTEG  RKS 

STORAG  IC(3)  ,NUM1),DEN<4) 

TABLE  IC(l-3)=4*0,M)M(l)=1.0,DEN(l-4)=l.,10.f2*-.,0. 

CONTRL  FINTIM=10.,DELT=.l 

PARAM  IN  =  1.  ,F=1.,K=50. 

*  SOLUTION  0*=  BASIC  SYSTEM  FCR  TIME  RESPONSE 

ERR0R1=IK-F*CUT 
0UT=TPNFp<0t3,IC,NUMfDENfK*FRRnRl) 

*  SOLUTION  OF  MODEL  POP  SENSITIVITY  COEFFICIENTS 

*  SENSK=SENSITIVITY  TO  FORWARD  GAIN  K 

*  SENSF=SENSITIVITY  TO  FEEDBACK  GAIN  F 

SFNSK=CLT+SFNSF 

X=TRNFR(0»3,IC,NUM,DEN,K*SCNSK) 

SENSF=X*(-F) 

PRINT  .1.0UT,SENSK,SFNSF 

PREPAR  .l.OL'T.SFNSK.SENSF 

GRAPH  SAME ,TTMF, CUT, SENS K,SENSF 

PRPLOT  ONLY 

ENH 

STOP 
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APPENDIX  P 


DSL/36C  DIGITAL  SIMULATION  FRCGRAM  FOR  MULTILOOP  SYSTEM 
TITLE  MULTILCCP  SFRVO  SENSITIVITY  COEFFICIENTS 


INTEG  RKS 

CONTRL  FINTIV=1C.0,DELT=.01 

PARA*  Kl=40.  ,K2=1.  ,K3  =  1.  ,KA  =  2. ,FGAIN=1.,IN  =  1 
TAU4=1. 


,TAUl=4.t 


*  SOLLTICN  CF  VLLTILCCP  SYSTFV  FCR  TIME  RESPONSE 

ERRCR1=IN-PGAIN*CUT 

D=EXCI7*(-K4)*TAU4 

FED1=D 

FED2  =  R«=ALPL<0. ,TAU4,0) 

FED*=FEC1-FEC2 

FPROR?=FPRCRl+FEDB 

EXCIT  =  RFALPL(0.  ,TAU1 fKl*ERR0R2) 

GFN=PCALPL<0.,2. ,K2*EXCIT) 

OUT=RFALPL<0. ,1. ,K3*GFN) 

*  SOLUTION  CF  MCOEL  FOR  SENSITIVITY  COEFFICIENTS 


SFN 

nx= 

FED 
FED 
<:PN 

SrN 
SFN 
P=- 
Q=R 
SFN 
X=R 
Y=R 
Z  =  R 


SK3=CLT-FGAIN*7 

X*<-K4)*TAU4 

XI  =C  X 

X2=PFALPl.(0.  ,TAU4,CX) 

SK4=FEDX1-FEDX2 

ST4=REALPL<0. ,TAU^  ,SCNSK4 ) 

SK1=SENSK3+SENSK4 

SPNSK! 

FALPL(0.,TAU1  ,P) 

ST1=P-Q 

FALPL(f).  ,4.  ,K1*SENS*1) 

EALPL(C.,2.,K2*X) 

»=ALPL(0.  tl«  ,K3*Y) 


PPINT   .1,0UT,SENSK3,SENSK1  ,SFNSU  f  SFNSK4,  SENST4 

PREPAR  .l,OUT,SENSK3.SFNSKl .SPNST11SENSK4.SENST4 

GRAPH  SAME,  TIME, CUT. SENSK3,SENSKlfSENSTl , SENSK4 , SFNST4 

GRAPH  SAME,TIME,SENSK1 tSENSTl 

PRPLOT  ONLY 

END 

STOP 
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APPFNDIX  C 


DSL/36C  DIGITAL  SIMULATION  PROGRAM  FOR  SIXTH  ORDER 
DIFFERENTIAL  EQUATION 

TITLE  SIXTH  ORDER  DIFFERENTIAL  EQUATION  SENSITIVITY  MODEL 

INTEG  RKS 

CONTRL  FINTIM  =  50..DELT  =  .Oe5 

PAR  AM  A  =  l.  ,B  =  3.4  ,  C  =11.4,0=9.  4  fE=7. 4  ^  =  .575,0*. 452 

*  SOLUTION  OF  DE  FCR  TIME  RESPONSE 

IN=STFF(C) 

XD6  =  (-B*XD5H-(-C*XD4)  +  (-D*XD3)-H-E*XD2)  +  (-F*X01)... 

♦(-G*X)+IN 

XD5=INTGRL(0.,XD6) 

XD4=INTGRL(0.,XD5) 

XD3=INTGRL(0.,XD4) 

XD2=IN7GPL<0. ,XD3) 

XD1=INTGPL(0.,XD2) 

X=INTGRL(0. ,XD1) 

*  SOLUTION  OF  MODEL  FOR  SENSITIVITY  CCFFFICIENTS 

SENSG=-G*SX 

SENSF=-F*SXD1 

SENSF=-F*SX02 

SENSD=-D*SXD3 

SCNSC=-C*SX04 

SXD6=SENSB*SENSC  +  SENSD  +  SFNSF+SEKSF-i-SENSG+X 

SXD«5  =  INTGPL(0. ,SXD6) 

SXD^=INTGRL(C. ,SXD5) 

SXD3=INTGPL(C. ,SXD4) 

SXD2=IKTGRL(0. ,SXD3) 

SXD1=INTGRL(C. fSXD2) 

SX=INTGRL(C. ,SXD1» 


NSE, SENSE, SENSG 
NSE,SENSF, SENSG 
NSO,SENSE,SENSF,SENSG 
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APPENDIX  D 


DSL/360  DIGITAL  SIMULATION  PRCGRAM  FOR  NONLINEAR  STSTEM 

TITLF  NONLINEAR  SYSTEM  SENSITIVITY 

INTEG  RKS 

STORAG  IC(3),B(4I 

TABLE  IC(l-3)=3*0,B(l-4)=l.,lO.,24.,O. 

CONTRL  FINTIM=12. ,DELT=.l 

PARAM  K  =  50. ,P1=-. 2, P2  =  . 2 

*  SOLLTICN  CF  NCM.INFAP  SYSTEM  FOR  TIME  RESPONSE 

IN=STFP(C> 

ERROR=TN-CLT 

Y  =  LIMI7(P1  ,P2, ERROR) 

0UT=DECM3fICf8,K*YJ 

*  SOLUTION    CF    WOCEL    FOR    SENSITIVITY    COEFFICIENTS 

*  SENSK=SENSITIVITY    TO    FORWARD    GAIN    K 

NOSORT 

IF(EPRCR.GT.P2.CR.ERRCR.LT.P1)GC    TO    1 

M=OLT-N 

GO    TO    2 

1  m  =  C 

2  CONTINLE 
N=OEQN(3,IC,B,K*M) 
SENSK=CLT-N 

PRINT  ,1,FRRCP,CLT,SENSK 

PREPAR  .l,OLT,SENSK 

GRAPH  SAVf ,TIMEtCUTfSENSK 

PRPLOT  ONLY 

END 

STOP 
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